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Abstract 

First, an algebraic criterion for integrability is discussed -the so-called 'superintegrability'- 
and some results on the classification of superintegrable quantum spin Hamiltonians based 
on si (2) are obtained. 

Next, the massive phases of the Z n -chiral Potts quantum spin chain (a model that violates 
parity) are studied in detail. It is shown that the excitation spectrum of the massive 
high-temperature phase can be explained in terms of n — 1 fundamental quasiparticles. 
We compute correlation functions from a perturbative and numerical evaluation of the 
groundstate for the Z3-chain. In addition to an exponential decay we observe an oscillating 
contribution. The oscillation length seems to be related to the asymmetry of the dispersion 
relations. We show that this relation is exact at special values of the parameters for general 
Z n using a form factor expansion. 

Finally, we discuss several aspects of extended conformal algebras (W-algebras). We ob- 
serve an analogy between boundary conditions for Z n -spin chains and W-algebras and then 
turn to statements about the structure of W-algebras. In particular, we briefly summarize 
results on unifying structures present in the space of all quantum W-algebras. 
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Mathematik und Physik sind die beiden theoretischen Erkenntnisse der Vernunft, welche ihre Objecte 
a priori bestimmen sollen, die erstere ganz rein, die zweite wenigstens zum Theil rein, dann aber auch 
nach Massgabe anderer Erkenntnisquellen als der der Vernunft. 

Immanuel Kant, Vorrede zur zweiten Auflage der "Kritik der reinen Vernunft" 

1. Introduction 



Theoretical physicists have paid much attention to low dimensions in the last years. On the 
one hand, models in low dimensions can more easily be treated also on a mathematically 
rigorous level, but on the other hand, low dimensions have indeed turned out to be relevant 
in the description of important experiments. This thesis is centered around two aspects 
of low-dimensional theoretical physics: Z n -quantum spin models and extended conformal 
algebras. 

Among the Z n -quantum spin models we will mainly concentrate on a 'Z n -chiral Potts 
model'. The first chiral Potts model that was introduced in 1981 by Ostlund in order 
to describe incommensurate phases of physisorbed systems [114] was a classical 2D spin 
model. The associated quantum chain Hamiltonians were obtained in 1981-82 by Rit- 
tenberg et al. [108,39]. Because this chain was not self-dual the location of the critical 
manifold was difficult. In 1983, Howes, Kadanoff and den Nijs introduced a self-dual 
Z3-symmetric chiral quantum chain [97], which however, does not correspond to a two- 
dimensional model with positive Boltzmann weights. Soon afterwards, von Gehlen and 
Rittenberg noticed that the remarkable property of the first gap of this model being lin- 
ear in the inverse temperature also applies to the second gap and can be generalized to 
arbitrary Z n [70]. Furthermore, the authors of [70] showed that the Ising-like form of 
the eigenvalues is related to this Z n -Hamiltonian satisfying the Dolan-Grady integrabil- 
ity condition [45] which is equivalent [115,42,43] to Onsager's algebra [113]. It was then 
shown by Au-Yang, Baxter, McCoy, Perk et al. that this integrability property - nowadays 
called 'super integrability' - can be implemented in a 2D classical model with Boltzmann 
weights defined on higher genus Riemann surfaces that satisfy a generalized Yang-Baxter 
relation. In the sequel the chiral Potts model attracted much attention because of these 
mathematical aspects, i.e. on the one hand due to the generalized Yang-Baxter relations 
[4, 5, 6, 10, 11, 12, 21, 22, 109, 110, 119, 121] and on the other hand because of Onsager's alge- 
bra [42, 43, 120, 1]. We will present new results showing that the model is also 'physically' 
very interesting although it is not directly related to a realistic 2D physisorbed system. 

We start by a general discussion of quantum spin models and indicate how one might 
systematically construct integrable nearest neighbour interaction quantum spin models 
exclusively by algebraic means. Afterwards, we specialize to Z n -quantum spin chains but 
drop the requirement of integrability. 

One of our main interests will be the excitation spectrum of Z n -quantum spin chains. 
We will show that in the massive phases the excitation spectrum can be explained in 
terms of quasiparticles. Although this is not very surprising, a particle interpretation is 
not directly incorporated into the spin chain Hamiltonian. Therefore it was not recognized 
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until very recently that a quasiparticle interpretation of the complete excitation spectrum 
is possible. Using numerical techniques we were able to show that the low-lying excitations 
in the zero momentum sectors can be explained in terms of n — 1 fundamental particles 
for n = 3, 4 at general values of the parameters [68, 66] and checked for n = 3 that 
this quasiparticle picture extends to general momenta [69]. For the superintegrable Z 3 - 
chiral Potts model McCoy et al. have derived a quasiparticle picture of the complete 
spectrum using techniques related to the Bethe ansatz [41]. Recently, they argued that this 
quasiparticle picture should in general be valid for the integrable Z3-chiral Potts quantum 
chain [103]. In this thesis we will show that both results can be combined into the general 
statement that the massive high-temperature phases of general Z n spin quantum chains 
have quasiparticle spectra. In fact, this quasiparticle picture will in certain cases give small 
corrections to the additivity of energy in the momentum zero sectors observed in [68]. 

Another main issue will be correlation functions of the Z n -chiral Potts model. Since 
a scaling exponent for the wave vector in the low-temperature phase of the Z3-chiral Potts 
model has been calculated in [105, 66] from level crossings in the ground state one expects 
an oscillatory contribution to the correlation functions also in the massive phases. We 
will demonstrate the presence of this oscillation by computing low- and high-temperature 
series and numerical evaluation of the ground state. These computations are limited to 
short distances. Therefore we also use a form factor expansion for the correlation functions 
which leads to a better conceptual understanding of the oscillation. 

These Z n -spin quantum chains have conformally invariant second order phase transi- 
tions for particular values of the parameters. They can be described by rational models 
of Z n parafermions [56] with c > 1 (for n > 4). The Z n parafermion algebras have the 
disadvantage that they are non-local and that each algebra only describes precisely one 
rational conformal field theory. However, Virasoro-minimal models have c < 1 [24] and are 
therefore not useful for describing these second order phase transitions. The discovery of 
extended conformal algebras [129] (also called W-algebras - for a review see e.g. [33]) was a 
major simplification because these algebras are local and -like the Virasoro algebra- admit 
infinite series of minimal models. Nowadays the second order phase transitions of Z n -spin 
chains can conveniently be described by the first unitary minimal models of the so-called 
WL4 n -i-algebras [8,57,54,55]. Similarly, one also needs c > 1 for applications in string 
theory. Very important applications require N = 2 supersymmetric rational conformal 
fields theories with c = 9 (see e.g. [75, 76]). Consequently, the classification of all rational 
conformal field theories is a natural question. Since those with c > 1 are obtained from 
W-algebras one may expect that W-algebras play an important role in this classification 
program. Another motivation for the study of W-algebras is that W-algebras have a rich 
mathematical structure, their complete classification still being an open question. 

With this thesis we will also try to contribute to these issues. Once again, we will look at 
the aforementioned second order phase transitions from the point of view of representations 
of W-algebras with non-trivial boundary conditions. We then turn to questions relevant 
for the classification problem. The classical versions of W-algebras are (to some extent) 
more easy to handle. We will therefore use them to examine cosets and orbifolds after 
addressing some basic problems of covariance with respect to the Virasoro subalgebra of 
the classical W-algebra. Some puzzles (see e.g. [51]) and also some confusion concerning 
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the structure of quantum W-algebras have recently been explained with the concept of 
'unifying W-algebras' [27,28]. We summarize these results and discuss their relevance for 
the classification problem. 

Perturbations of conformal field theories (see e.g. [Ill]) describe the off-critical scaling 
region in statistical mechanics. For the Z n -chiral Potts models two types of perturbations 
are relevant: A thermal perturbation (see e.g. [130, 131, 132, 111]) that makes the models 
massive and a chiral perturbation [38] that breaks parity. Both perturbations separately are 
integrable but it is still an open question if the combined perturbation is also integrable. 
W-algebras might also be useful for addressing problems of this type which is however 
beyond the scope of this thesis. 

The outline of this thesis is as follows. Chapter 2 contains a review of the algebraic 
formulation of superintegrability and indicates how superintegrable quantum spin models 
can be constructed systematically. In chapter 3 we summarize basic facts about the chiral 
Potts model and perturbation series. In particular, in section 3.2 we present a precise 
duality statement. Then, the excitation spectrum of Z n -spin chains is examined in chapter 
4. After reviewing the quasiparticle picture for the superintegrable Z3-chiral Potts model 
[41] we start by studying the low-lying levels with high- and low-temperature series. These 
results illustrate a general quasiparticle picture which we discuss in section 4.4. This 
quasiparticle picture is strictly valid only in the infinite lattice limit but we also discuss 
finite-size corrections. Since we use perturbative arguments for the derivation of this 
quasiparticle picture we also have to discuss the radius of convergence in section 4.5. The 
various aspects of the correlation functions mentioned above are discussed in chapter 5. 

The final chapter 6 shifts to a different aspect of low-dimensional theoretical physics: 
Extended conformal algebras. In section 6.2 we show how automorphisms of Casimir W- 
algebras can be used to find extensions of rational conformal field theories by imposing 
twisted boundary conditions. The relevance of these results to Z n -quantum spin chains is 
discussed in section 6.3. We then turn to classical W-algebras in section 6.4 and discuss 
their stt(l, 1) structure. For classical W-algebras one can get a good handle on orbifolds 
and cosets as we show in section 6.5. These results are useful also for obtaining new 
insights in the structure of quantum W-algebras. Some new statements about the space 
of quantum W-algebras are summarized in section 6.6 where in particular the notion of 
'unifying W-algebras' is explained. 
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This thesis is (to a large extent) a summary of already published material. More precisely, 
the following sections are based on the following publications: 



Sections 3.2., 4.3. and 5.3.: 

N.S. Han, A. Honecker, Low- Temperature Expansions and Correlation Functions of 
the Z 3 -Chiral Potts Model, J. Phys. A: Math. Gen. 27 (1994) p. 9 

Section 4.2.: 

G. von Gehlen, A. Honecker, Multi-Particle Structure in the Z* n -Chiral Potts Models, 
J. Phys. A: Math. Gen. 26 (1993) p. 1275 

Sections 4.2., 4.4., 4.5., 5.1., 5.4. and appendices A, B: 

A. Honecker, A Perturbative Approach to the Chiral Potts Model, preprint BONN- 
TH-94-21, frep-th/940912l 



Section 5.2.: 

G. von Gehlen, A. Honecker, Excitation Spectrum and Correlation Functions of the 
Tj^-Chiral Potts Quantum Spin Chain, Nucl. Phys. B435 (1995) p. 505 

Sections 6.2. and 6.3.: 

A. Honecker, Automorphisms of W '-Algebras and Extended Rational Conformal Field 
Theories, Nucl. Phys. B400 (1993) p. 574 

Section 6.5.: 

J. de Boer, A. Honecker, L. Feher, A Class of W- Algebras with Infinitely Generated 
Classical Limit, Nucl. Phys. B420 (1994) p. 409 

Section 6.6.: 

R. Blumenhagen, W. Eholzer, A. Honecker, K. Hornfeck, R. Hiibel, Unifying W- 
Algebras, Phys. Lett. B332 (1994) p. 51 

R. Blumenhagen, W. Eholzer, A. Honecker, K. Hornfeck, R. Hiibel, Coset Realization 
of Unifying W- Algebras, preprint BONN-TH-94-11, DFTT-25/94, [hep-th/ 9406203] , to 
appear in Int. Jour, of Mod. Phys. A 

There are some minor extensions throughout these sections. The previously unpublished 
material is concentrated mainly in sections 2.2. - 2.4. and section 6.4. 
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/ was at the mathematical school, where the master taught his pupils after a method scarce imaginable 
to us in Europe. The proposition and demonstration were fairly written on a thin wafer, with ink 
composed of a cephalic tincture. This the student was to swallow upon a fasting stomach, and for 
three days following eat nothing but bread and water. As the wafer digested the tincture mounted to 
the brain, bearing the proposition along with it. 

Jonathan Swift in "Gulliver's Travels" 

2. Superintegrable quantum spin models 

By now an abundance of integrable classical spin models in two space dimensions, or 
equivalently quantum spin chains in one space dimension, are known. The quantum spin 
chain Hamiltonians in turn can be regarded as evolution operators of a lattice field theory in 
one plus one dimensions. This indicates that it would be desirable to know also integrable 
spin models in higher dimensions in order not to be restricted to surface physics or other 
planar media. Unfortunately, very little is known so far about higher dimensions despite 
major efforts. 

A promising idea was put forward in [45] using some algebraic criterion for a particular 
kind of integrability - the so-called ' super integr ability'. To our knowledge, no results 
in dimensions greater than one have been obtained with this criterion yet either. In 
this section we would like to give some indications how one can systematically construct 
superintegrable quantum spin Hamiltonians. So far, we have not yet reached the goal of 
finding integrable models in higher dimensions. Nevertheless we show that the problem 
is in principle amenable to computer computation and we also obtain a few related new 
results 1 ). 

2.1. The general setting 

Def. I: (Nearest neighbour interaction quantum spin model) Let A be a lattice of dimension d 
and T be its fundamental domain. Furthermore, consider an associative algebra of dimension 
n with generators fi* including the identity 1. For any finite subset S(A) C A with p points 
define operators fit associated to 'site' x e S(A) by operators l®...®l®fi l (g>Ilcg>...Cg>l 
in the p-fold tensor product of the algebra with itself. Then, a general nearest neighbour 
interaction quantum spin Hamiltonian in d dimensions is given by 

n n n 

^(general) = £ £ ^ + A £ £ £ £ vfjQ^. (2.1.1) 
x€S(A) i=l xES(A) eEJ 7 i=l 3 = 1 

A Hamiltonian of this type is called 'isotropic' if the coupling constants vf do not depend 
on the direction e of the interaction, i.e. vf „■ = u» , for all e. 

The Hamiltonian (2.1.1) is of the form 

H = H + XV. (2.1.2) 



x ) Note that recently a different algebraic approach to the construction of integrable 
quantum spin models in one dimension has been proposed in [81]. 
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It was shown in [45,42,43] that Hamiltonians of this form are integrable if 

[H , [H , [Ho, V]]] = Cl [H , V] , (2.1.3a) 
[V,[V,[V,Ho]]]=c 2 [V,H ]. (2.1.3b) 

We will call a Hamiltonian of type (2.1.2) satisfying (2.1.3) ' superintegrable^ (note that 
this terminology is not used entirely consistent in the literature - in contrast to us, some 
authors include the generalized Yang-Baxter relations in the notion of superintegrability). 

Clearly, one can modify c\ and c 2 by rescaling the individual terms Ho and V of the 
Hamiltonian. This is usually exploited in order to set c\ = c 2 = 16. However, we will not 
impose this normalization condition but will instead normalize Hq and V suitably. One 
should keep in mind that if Ho + XV is superintegrable, so is (hHo) + X(vV) for arbitrary 
h, v. 

The condition (2.1.3) is trivially satisfied if [Hq, V] = 0. In this case, one can simultane- 
ously diagonalize Hq and V and in this manner diagonalize H for general A. Clearly, for 
[Hq, V] = all eigenvalues of H are linear in A. Therefore, if A should play the role of tem- 
perature, the Hamiltonian H = Ho + XV gives rise to trivial thermodynamics. Thus, the 
non-trivial physics must be given by the potential, usually for A — > oo, i.e. the physically 
relevant Hamiltonian is H = V + XHq with small A. In this case the diagonalization of V 
may still be a difficult task although the conservation of Hq can be very helpful. 

For [Hq, V] 7^ one can recursively generate Onsager's algebra [113] from A = Hq and 
A 1 = V (see [45,42,43]): 

[A m -, A n ] = 4G m — n , [G m , A n ] = 2A n j rm — 2A n _ m , [G m , G n ] = (2.1.4) 

if (2.1.3) is satisfied with c\ = c 2 = 16. From (2.1.4) one obtains the following infinite set 
of conserved commuting charges: 

-Qm = \ (A m + A_ m + XA m+1 + XA_ m+1 ) . (2.1.5) 

Note that G_i = -G x = \[Hq,V\. Thus, if [H ,V] = the recurrence relations (2.1.4) 
can be trivially solved leading to A m = A n if m = n mod 2. In this case, one has Q m = Q n 
for m = n mod 2, i.e. the Hamiltonian is not necessarily integrable in the sense that it has 
an infinite set of conserved commuting charges. This motivates the following definitions: 

Def. II: A Hamiltonian of the form (2.1.2) is called 'quasi superintegrable if it satisfies 

[H ,V]=0. 

Def. Ill: A Hamiltonian satisfying (2.1.3) with [Hq, V] ^ is called 'strictly superintegrable . 
2.2. First considerations 



We will be interested in associative algebras that have finite dimensional (complex) defining 
representations. These algebras can be considered as subalgebras of gl(n, C) for some n. In 
the Hamiltonian (2.1.1) we may restrict to operators fi l G sl(n, C) for the following reason: 
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The complement of sl(n, C) in gl(n, C) is spanned by the identity 1. If the identity 11 should 
appear in the potential term V, this part of the potential could be absorbed in the free 
part Hq. Now, if the identity would appear in the free part Hq it would shift the ground 
state energy per site eo by a constant amount and not change the excitation spectrum. 
Such a modification may be neglected because of its physical irrelevance. Thus, we assume 
Q i G sl(n,C) in (2.1.1). 

Next we note that the conditions (2.1.3) are local ones. We may therefore choose a suitable 
domain S(A) and suitable boundary conditions. Furthermore note that the eigenvalues of 
the Hamiltonian (2.1.1) will not depend on the details of the lattice A but only on the 
way how neighbours are linked. We will therefore restrict to cubic lattices in d dimensions. 
This leads us to the following specialization of the Hamiltonian (2.1.1): 

n 2 -l d n 2 -ln 2 -l 

^(general') = £ £ d^*, + A £ £ £ £ (2-2.1) 

seZ d i=1 seZ, d i=1 a=1 6=1 

where n is now the dimension of the defining representation of sl(n, C). On (2.2.1) we will 
impose periodic boundary conditions, i.e. the Xi will be identified if they are equal mod 
Ni. 

Finally we make some observations that will be useful for treating the problem on the 
computer. For a Hamiltonian of type (2.1.1) the conditions (2.1.3) will involve at most 
four-fold products of generators Qt sitting on four neighbouring sites. Due to the locality 
of the conditions (2.1.3) this means that we may choose all iVj = 4 when checking them 
for a general ansatz of the form (2.2.1). 

Note that the condition [Hq, V] = does not involve the dimensionality of the lattice. 
Therefore, for any solution {di, f a ,&} with [Hq, V] = in d = 1 we immediately obtain the 
corresponding general solution in arbitrary dimensions d by fixing the parameters {di, v J a b } 
such that for fixed j they satisfy the equations for d = 1 independently. 

After these first general considerations the problem of solving (2.1.3) is now straightforward 
although not necessarily very easy. 

2.3. The quasi superintegrable case 

As we have pointed out before we are mainly interested in strictly superintegrable spin 
chain Hamiltonians. However, in order to single them out among all solutions to (2.1.3) 
we have to find the quasi superintegrable ones first. 

Fortunately, this problem is comparably simple. In order to evaluate [H , V] for a Hamil- 
tonian of the form (2.2.1) we do in fact not need that the Q a form an associative algebra, 
we need nothing else than the Lie bracket. Introducing structure constants by 

[n a ,n b ] = J2fc b n c (2.3.1) 
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one obtains 

d n 2 -l 

[ H ^y]= E E E d «HJd +vhfc b }^Phe i - (2.3.2) 

Se Z d i=1 a,b,c,d=l 

Clearly, H + XV will be quasi superintegrable iff 

n 2 -l 

E d * Wjd + <6/c b } =0 Vc, d, i . (2.3.3) 

a, 6=1 

The condition (2.3.3) is very reminiscent of an invariance condition with respect to the 
generators of the Lie algebra for which d a ^ 0. Thus, it will not be surprising if the quasi 
superintegrable spin chain Hamiltonians turn out to have extra continuous symmetries. 



2.4. Results 



In this section we will present some new results for quantum spin models based on si (2, C) 
and afterwards recall some facts about Z n spin chains. After complexification s/(2,C) is 
the same as su(2). We can therefore use the standard Pauli spin matrices 

i -°o. ■■-(! ;)• -(! o)- < 2 -» 

In fact, we will not need the matrix realization given by (2.4.1). Instead we will look at 
the associative algebra generated by 1, cr z , cr x , cr y , e.g. (<j z ) 2 = (cr x ) 2 = (a y ) 2 = E. 

In the basis (2.4.1) the isotropic Ising quantum spin model on a cubic lattice in d dimensions 
which is one special case of (2.2.1) is given by: 

d 

Rising) = _ ^ ff z_ x J- E^+eV (2A2) 

Note that superintegrability is known since Onsager's days for _£/"( Ism s) only in d = 1. For 
d = 2 we will show below that _£/"( Ism g) is not superintegrable. 

Let us now turn to the most general ansatz for a nearest neighbour interaction Hamiltonian 
based on sl(2, C) in d dimensions. In addition to the previous comments we recall that 
SO (3) acts on su{2) by inner conjugations. This can be exploited in order to fix the basis in 
sl(2, C) and the normalization of H such that H = — J2x a s- Thus, the general spin-1/2 
Hamiltonian becomes 

d 2 2 

H = - E ^-a E EEE<M + e ; (2-4.3) 

* e Z d a e Z d *=1 «=0 6=0 
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with 0° = a z , 1 = a x and O 2 = a y . Note that there may be some residual freedom of 
choice of basis corresponding to rotations around the z-axis. 

First we consider one-dimensional models. In d = 1 the general ansatz (2.4.3) becomes 

N N 2 2 

^ = "E^- A EEE (2-4.4) 

x=l x=l a=0 6=0 

With this ansatz one can show by a tedious but straightforward calculation that there are 
only two different inequivalent families of solutions to eq. (2.1.3) given by theorems I and 
II. 

Theorem I: The most general quasi superintegrable Hamiltonian of type (2.4-4) in d = 1 
is given by: 

N N 

Hil) = -J2<- X H W+i + + « W+i - + ™>x+i } ■ (2- 4 -5) 

x=l x=l 

m 

Remarks: 

1) The Hamiltonian (2.4.5) is invariant with respect to rotations in the two-dimensional 
space spanned by a x , a y . For v = v = this solution contains the Hamiltonian 
generated by [1] with (in the notations of [1]) g\ = g-±. Note that the Hamiltonian 
eq. (11) of [1] is not superintegrable for general g± ^ g-\. 

2) In [79] a spin-1 Hamiltonian of the type H = J2 X $x • Sx+i was studied numerically 
for up to 22 sites. This Hamiltonian is a special case of (2.4.5) and demonstrates the 
usefulness of conservation of Hq for numerical computations. 

The second solution to (2.1.3) is strictly superintegrable. For the second solution in d = 1 
we have some freedom of rotation between a x and a y . After fixing this freedom we arrive 
at 

Theorem II: The most general strictly superintegrable Hamiltonian of type (2.4-4) ^ n 
d = 1 is given by: 

N N 

H{2) = ~ E < - X E K^+i + " W+i - } ■ (2-4.6) 

x=l x=l 

m 

Remark: The solution (2.4.6) contains the Ising model (2.4.2) in d = 1 for v = 0. 

The conditions for a general Hamiltonian (2.4.3) to be superintegrable are already very 
complicated for spin 1/2 in d = 1. Therefore, it does not make sense to start with a general 
ansatz in d > 1 but one has to restrict oneself to special cases. 

A trivial (i-dimensional strictly superintegrable Hamiltonian can be obtained from (2.4.6) 
by setting all v * b = for % ^ k with some fixed k. For v\ b one then simply takes 
the solution (2.4.6). However, this solution is trivial since it consists of decoupled one 
dimensional models. Therefore, we will not consider solutions of this type any more. 
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First, we note that the Hamiltonian _H"( Ism s) given by eq. (2.4.2) is not superintegrable in 
d = 2. Even more strongly, one can show that 

Theorem III: There is no strictly superintegrable Z 2 charge conserving Hamiltonian of 
type (2.4-3) in d = 2. Here Z 2 charge conservation means that v q 1 = Vq 2 = v\ = v\ = 
fori = 1,2. m 

However, from the preceding discussions we obtain the 
Corollary: The Hamiltonian 

# (M) = - E E EM^ + *+^ + *) 

seZ d s e Z d *=1 (2 4 7) 

is quasi superintegrable for general d. 

Current knowledge about strictly superintegrable spin quantum chains based on sl(n) with 
n > 2 is restricted to the result of von Gehlen and Rittenberg who constructed the so-called 
'superintegrable Z, n -chiral Potts model [70]. More precisely, one has 

Theorem IV: ([70]) The Z n -chiral Potts spin quantum chain defined by 

<» = - f (2H, + A g (l - icot I* r-^ (2.4.8) 

is strictly superintegrable. H and V are operators in the fundamental representation of 
sl(n). H is the Cartan generator of the principal sl(2), i.e. H = diag ( I ^ L , ■ ■ ■ , — ■ 
Furthermore, one has Y = Q^ Q E a ) + E- 1 where the sum runs over the n — 1 simple 
positive roots and 7 is the highest root. Ep is the root vector associated to a root [3 in a 
normalization such that E a for a simple root a has a single non-zero entry equal to one 
on the first upper off-diagonal. ■ 

Ahn and Shigemoto [1] realized that one can add terms of type \gi (A m+ i + A_ m+ i) with 
arbitrary parameters g\ to the conserved charges Q m in (2.1.5), still ensuring the commu- 
tativity of these charges. Furthermore, A_i = V — ^[H ,V] is obviously still a nearest 
neighbour interaction for Hamiltonians of type (2.1.1). Thus, adding |A (A m -i + A_ m _i) 
to (2.1.5) they arrived at the following corollary to theorem IV: 

Corollary: ([1]) The Z n -spin quantum chain defined by 



N 

< = -£ 

i=i 



2H, + £ (l -icot f) (ATfl-tf + A=*=j + -f) 
fc=l v J 



(2.4.9) 



is integrable for general A, A, i.e. it gives rise to an infinite set of commuting charges. The 
additional operator H is defined by H = (^a^) — E- 1 . 

Remark: The Hamiltonian (2.4.9) is neither quasi nor strictly superintegrable for general 
values of the parameters. 
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Einstweilen wissen wir noch gar nicht, in welcher Sprache wir uber das Geschehen im Atom reden 
konnen. Wir haben zwar eine mathematische Sprache, das heiBt ein mathematisches Schema, mit Hilfe 
dessen wir die stationaren Zustande des Atoms oder Ubergangswahrscheinlichkeiten von einem Zustand 
zu einem anderen ausrechnen konnen. Aber wir wissen noch nicht — wenigstens noch nicht allgemein — 
wie diese Sprache mit der gewdhnlichen Sprache zusammenhangt. 

Werner Heisenberg in "Der Teil und das Ganze", Kapitel zu 1925-1926 



3. The chiral Potts quantum chain 



This chapter and chapters 4 and 5 of this thesis focus on the 'Z n -chiral Potts quantum 
chain' which is a generalization of (2.4.8). From here on we will rely on methods that do 
not make use of any particular integrability properties. 



3.1. Preliminaries 



This section summarizes well-known basic facts about Z n -spin quantum chains. We also 
introduce some notions that will be useful later on. For more details see e.g. the review 
[40]. 

In the following we will study the Z n -spin quantum chain with N sites which is defined by 
the Hamiltonian: 

N n-1 

h n } = - E E + x ^ r J r l+i} ■ ( 3 - L1 ) 

3 = 1 k=l 

For the low-temperature phase (large A) it is more convenient to use a different normal- 
ization than (3.1.1): 

N n-1 

= \hP = - e E + a ^ T 7+i} ( 3 - L2 ) 

3=1 k=l 

with A = A -1 . 

Oj and Tj freely generate a finite dimensional associative algebra with involution by the 
following relations (1 < j,l < N): 

o 3 o x = cncTj , r^r, = TtTj , ajTt = r^cA* , a? = = 1 , 
where u is the primitive nth root of unity u = e 2 ^ 1 . 

For the Hamiltonian (3.1.1) we also have to specify boundary conditions, i.e. we have to 
define IV+i- The following choices will be considered (in the terminology we follow Cardy 
[37]): 

1) Free boundary conditions 

riv+i = 0. (3.1.4a) 
11 



2) Toroidal boundary conditions, i.e. we identify the N + 1st site with the 1st site. 

a) Cyclic boundary conditions 

T N+1 = u- R T 1 . (3.1.4b) 

b) Periodic boundary conditions which are the R = special case of cyclic boundary 
conditions 

r N+ i = r 1 . (3.i.4c) 

c) Twisted boundary conditions 

T n+1 =uj- r T+. (3.1.4d) 
We will mainly focus on periodic boundary conditions. 

The Hamiltonian (3.1.1) contains 2n — 1 parameters: The temperature-like parameter A 
that we choose to be real and the complex constants ak and a^. is hermitean iff 

«fc = a* n _ k and ctk = ot* n _ k . 

The algebra (3.1.3) is conveniently represented in 

H N := C n ®C n ® ...® C n y (3.1.5) 

N times 

labeling the standard basis of C n by {e , . . . , e n _i}. Then a basis for (3.1.5) is given by: 

| i x . . . in) '■= <S> • • • <S> e iN , < ij < n - 1. (3.1.6) 

Now the following operation in the space (3.1.5) is a faithful irreducible representation r 
of the algebra (3.1.3): 

r(<Tj) \ i 1 ...i j ...i N ) = uj 1 * \h...ij...i N ) , ^ i ^ 

r(Tj) | i\ . . . ij . . .in) =\ii ■ ■ ■ (ij + 1 mod n) . . An) ■ 

The involution is the adjoint operation with respect to the standard scalar product in the 
tensor product of C n . For low-temperature expansions of (3.1.2) it is more convenient to 
consider a representation f different from r: 

r(Tj) \i 1 ...i j ...i N ) =u l > \h...ij ...i N ) ^ 1 g ^ 

r(o~j) | i\ . . . ij . . . %n) =| ii . . . (ij — 1 mod n) . . . ijv) . 

The Hamiltonian (3.1.1) commutes with the Z n charge operator Q := YljLi a j acting 
on the vectors (3.1.6) as 

r(Q) \h...i N ) =a;(^r=i i \i 1 ...i N ) (3.1.9) 

which shows that the eigenvalues of Q have the form iv® with Q integer. Thus, has 
n charge sectors which we shall refer to by Q = 0, . . . , n — 1. 
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also commutes with the translation operator T N that acts on the basis vectors 
(3.1.6) in the following way: 

r(T N ) | ziz 2 . . . i N ) =|*2 • • • iNh) ■ (3.1.10) 

The eigenvalues of T/v are A^th roots of unity. We label them by e lP and call P the 
'momentum'. We choose < P < 2n corresponding to the first Brillouin zone and have 
PG{0,|,..., 27t( n~ 1} }- Note that the states 

\\iii 2 . . .i N -ii N )) P := -)= ( \ hi 2 . . .ijv-iijv) + e tP | Wi*2 ■ • -ijv-i) + • • • 

^ V (3.1.11) 
+ e^" 1 ) |i 2 ...x JV -i<i V *i>J 

are eigenstates of T/v with eigenvalue e tP '. A/" is a suitable normalization constant. If the 
state I zi . . . in) has no symmetry (i.e. T N \ ii . . .i N ) 7^ ^ . . . i N ) for all < k < N), one 
has M = N. This will apply to most cases below where we need (3.1.11). For a definition 
of the momentum eigenstates (3.1.11) in case of cyclic boundary conditions R 7^ see e.g. 
[73,40] 

We will frequently use the following parameterization of the constants aj~ and a^, 
fixing their dependence on k: 

"A; = — . nk , OL k = — , . • (3.1.12) 
sin — sin — 

n n 

This is a suitable choice because it includes a large class of interesting models. 

For = (p = one obtains real a k = &k = ■ htn ■ This leads to models with a second order 

phase transition at A = 1 which can be described by a parafermionic conformal field theory 
in the limit iV — > 00 at criticality [56,8]. These so-called Fateev-Zamolodchikov-models 
[57] lead to extended conformal algebras WA n -i where the simple fields have conformal 
dimension 2, . . . , n [54]. The spectrum of the Hamiltonian (3.1.1) can be described by the 
first unitary minimal model of the algebra WA n -i- For n = 3 the symmetry algebra is 
Zamolodchikov's well-known spin-three extended conformal algebra [129] at c = |. 

Choosing <fi = ip = ^ in (3.1.12) for the Hamiltonian (3.1.1) one recovers the superinte- 
grable Z n -chiral Potts model (2.4.8). 

The parameterization (3.1.12) also includes the family of integrable models discovered in 
refs. [4, 5, 6, 10, 11, 12] which interpolates between the integrable cases at <p = = 0, A = 1 
and (j) = <p = The Hamiltonian (3.1.1) is integrable if one imposes the additional 
constraint 

cos if = A cos (p (3.1.13) 

(or equivalently A cos <p = cos cj) for the low-temperature phase) on the parameterization 
(3.1.12). For <j) = if = this yields A = 1 - the conformally invariant critical points. At 
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= ip = 0, the Hamiltonian is self-dual, i.e. it is invariant under a duality-transformation 
such that Hj^\\) = XH < ^\\~ 1 ). The Hamiltonian is also self-dual on the superintegrable 
line (f) = (p=z. Hff> with the choices (3.1.12), (3.1.13) is in general not self-dual any more 
whereas particular choices yield a self-dual Hamiltonian. If we choose for (3.1.12) = (p 

and neglect (3.1.13) will be self-dual again. Therefore we choose to refer to (3.1.1) 

with (3.1.12) as the general 'chiral Potts model'. We will not consider the integrable case 
where the additional constraint (3.1.13) is satisfied in detail. 




The main features of the phase diagram of the Z n -ch^ai P^ptts chain are already present 
for n = 3. The Z 3 version of (3.1.1) is known to have four phases [4,6,66]: Two massive 
and two massless phases (see Fig. 1). One of the massive phases is ordered and the other 
massive phase is disordered. In the following we will mainly be interested in the massive 
phases. At (p = f = ^ the massive high-temperature phase (small A) appears in the 
range < A < 0.901292 . . . and the massive low-temperature phase appears in the range 
1/0.901292 . . . < A < oo [4, 6]. There is a second order phase transition at the point A = 1 
for (j) = ip = and general Z n . This conformally invariant point is described by the models 
of Fateev and Zamolodchikov with WA n -i -symmetry [56,57,8,54]. The conformal field 
theories can be perturbed by a thermal operator (see e.g. [58,130,131,132,111]). This 
perturbation is integrable and corresponds to going to A ^ 1 on the line </> = </? = of 
the Z n -quantum spin chain. One can also introduce an integrable chiral perturbation of 
the conformal field theories [38]. The Z n -spin quantum chains describe both this chiral 
perturbation as well as the integrable thermal perturbations. 



14 



Our main interest is the spectrum in the limit N — > oo of Of course, we have 

to specify how the limit is to be taken. One possible definition was given in [91]. Here, we 
want to omit the technical details. We just mention that the basic idea is first to identify 
momentum eigenstates for chains of different lengths and then to consider the weak limit 
of T/v and the excitation spectrum operator 

AH^ := - E%n (3.1.14) 
where is the groundstate energy for N sites. 

3.2. Duality of spectra 

There is a first statement about the spectra of the Hamiltonian (3.1.1) which allows us 
to restrict to the high-temperature phase A < 1 for the discussion of the spectrum: The 
spectra in the high-temperature phase at A are dual to those in the low-temperature phase 
at A = A if we interchange ctfc and ctk- This statement for Z n -quantum spin chains has been 
known for a long time [52, 92] and was also used in [97]. However, special attention has to 
be paid to the boundary conditions when performing duality transformations. It has been 
observed in [70] that the duality transformation interchanges the role of the charge Q and 
boundary conditions R. The precise statement is given by the following 

Theorem V: (Duality) Denote the Hamiltonian (3.1.1) by H^(\, R ht , , a^) and 

(3.1.2) by H^\\, R lt , a l £ , a%) , including in both cases explicitly the corresponding param- 
eters. Furthermore, abbreviate the space with charge Q ht in the high-temperature phase by 
H9 M and the eigenspace of r(Q) to eigenvalues^ by W9 H . Then H^\\, R ht , a'jf , ajjj*) 
restricted to HP and H^\\ R lt , a l £, a l £) restricted to JiP have the same spectra if 

Q lt = R ht , R lt = Q ht , 

(3 2 Tl 

^lt _ n ht JLt _ frht \-\ V ' ' 

CHu — QLu , Q>u — CXi, , A — A. 



Proof: We present a simple non-standard proof of theorem V [82] that is different from 
the approach of e.g. [52,92]. We derive duality by comparing the representation r (3.1.7) 
to the representation f (3.1.8). Working with the states has the advantage that it is 
straightforward to take care of the boundary conditions. Note that r(Tj) = r(o~j) and 

3 ' 

state | GS) ^ to be the ground state in TiP 



r(o~j) = r(rj) and that the representations r and r are unitarily equivalent. Now fix the 

-.a 



n— 1 

GS) Q := —j= ^2 u l ' Q \ . (3.2.2) 



1=0 



Then the following states are a basis for HP U : 

\Q;i 2 -.-i N ) :=r(a- t2 )...r(a- tN ) |GS) Q ". (3.2.3) 
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Note that this implies: 

f (<xi) \ Q;i 2 ...i N ) = u QH | Q- (z 2 + 1) . . . (i N + 1)> . (3.2.4) 
Now consider the following intertwining isomorphism /: 

I \Q;i 2 . ■ -i N ) : = \(-i 2 )(i2 - 13) ■ ■ - (iN-i ~ iN)(iN + R H )) ■ (3.2.5) 
It is now straightforward to check using the basis (3.2.3) that 

Ir{au+i mod N) ) = r(r j r+ +1 )/ , Jf (r^) = r{a 3 )I . (3.2.6) 

The observation that / is a unitary map and r and f are unitarily equivalent in conjunction 
with (3.2.6) concludes the proof. □ 

Remarks: 

1) The momentum decomposition can be applied alike in the high- and low-temperature 
phase. Thus, theorem V is also valid if we further restrict to eigenspaces with mo- 
mentum P. 

2) Duality (3.2.1) preserves the condition (3.1.13). Thus, each integrable chiral Potts 
model is dual to exactly one integrable chiral Potts model. 

3.3. Generalities about perturbation theory 

Some of the results that will be reported later on rely on perturbation theory. Therefore, 
we now review the general outline for perturbation theory to all orders as presented in [23] 
which directly applies to the degenerate case as well. It should be clear to the reader that 
this method is not in all cases the most powerful one available. Nevertheless, unlike e.g. 
cluster expansions it can be directly applied to any problem we are interested in (including 
higher excitations and correlation functions). Note also that we do not intend to study 
the massless regimes III and IV of Fig. 1 for which perturbation expansions are certainly 
not well suited. 

The Hamiltonian (3.1.1) can be written as H = H + XV with H = — . k ®k&y V = 
-Ej,k a ^]+i- The eigenstates for H are obvious, thus we have solved: 

H I a) =E$ I a). (3.3.1) 

Now one can solve H | a(A)) = | a(A)) for small A as follows: Let be the 

projector onto the eigenspace of H with eigenvalue E^y We can treat non-degenerate 
and degenerate perturbation theory alike if we choose | a) such that 

?w v 1 °> = E S 1 a > ( 3 - 3 - 2 ) 

with a constant E^ , i.e. Vq^ is to be chosen diagonal. One also needs a regularized 
resolvent g(z) of H : 

g(z) := (TL-q^iz-Ho)- 1 . (3.3.3) 
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Then, the Wigner-Brillouin perturbation series 

oo oo 

E » = E yE t) ' i a ^ = E y i a ' ^ ( 3 - 3 - 4 ) 

are given by the following recurrence relations [23]: 
|a,0) =|a) 

| a, i/> =^° ) ) )|y \a,u-l) -J2\a^-f,)E$y (3.3.5) 

4> +1) = H K|a,i/>. 

Note that neither | a(A)) nor | a, z/) are in general normalized although | a) must be nor- 
malized to one. Observe that the derivation of (3.3.5) does not rely on H being hermitean. 
Therefore, (3.3.5) may also be applied to diagonalizable but non-hermitean H. 

The radius of convergence of the series (3.3.4) can be more easily discussed in a different 
framework. Therefore, we postpone such a discussion to section 4.5. 



There is one observation that makes explicit evaluation of high orders for the Z n - 
Hamiltonian (3.1.1) possible. The energy-eigenvalues E^ of Hffl do depend on the chain 

length N. However, for the low lying gaps A£| a } 2 ) of AHffl (see (3.1.14) ) the coefficients 
for powers of A become independent of iV up to order A^ -2 (see e.g. [7]). Intuitively, this 
can be inferred from the fact that (3.1.1) shows only nearest neighbour interaction and 
thus we need N — 1 powers in V to bring us around a chain of length N. Smaller powers 
in V (or A) do not feel that the length of the chain is finite. 



2 ) This will apply precisely to the fundamental quasiparticle states to be discussed below. 
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Paa den anden Side m0der vi Vanskeligheder af en saa dybtliggende Natur, at vi ikke aner Vejen til 
deres L0sning; efter min personlige Opfattelse er disse Vanskeligheder af en saadan Art, at de naeppe 
lader os haabe indenfor Atomernes Verden at gennemf0re en Beskrivelse i Rum og Tid af den Art, der 
svarer til vore Saedvanlige Sansebilleder. 

Niels Bohr, Letter to H0ffding in 1922 

4. Excitation spectrum 

In this chapter we will study the excitation spectrum of the Z n -chiral Potts model. For 
this end we denote the kth gap in charge sector Q by AEq^ starting to count with k = 0. 
More precisely, we define 

AE Q , k (P) := E Q , k (P) - £ ,o(0) (4.1) 

where Eq^{P) is the fcth energy level with momentum P in charge sector Q. In the massive 
phases the groundstate has Q = P = 0, i.e. i?o,o(0) is the energy of the groundstate. We 
will mainly be interested in the thermodynamic limit, i.e. in the limit N — > oo of (4.1). 
Thus, if we write down energy gaps AEQ jk (P) we will usually think of sufficiently large 
lattices. It will be clear from the context when gaps are evaluated for particular finite 
chain length N . 

In general we have the following symmetry under charge conjugation 

A£ n _ Qifc (P, </>, ip) = AE Q , k (P, -<f>, -<p). (4.2) 

Therefore, it is sufficient to restrict to those charge sectors with Q < ^ in all explicit 
calculations where we keep the dependence on and cp. 

4.1. Analytic results for the superintegrable Z 3 -chain 

The results in [41, 103] and [68, 69, 91] suggest that we may expect a quasiparticle spectrum 
in the complete massive high-temperature phase of the general Z n -Hamiltonian (3.1.1). At 
(f) = <p = ^ the spectrum of the Hamiltonian (3.1.1) has been determined analytically by 
McCoy et al. [41] for n = 3. Their result is that the spectrum is a perfect quasiparticle 
spectrum. We formulate their precise result in this section. 

Def. IV: ([41]) Let 



Ei(P r ) = e r = 2 | 1 - A 



E 2 ( J P 2s )=e 2s =4|l-A 



(4.1.1) 
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with the auxiliary variables v r and V2 S related to the momenta by: 

e iP r _ i 
u)(u)e lPr — 1) 



iP r _ -y 

v r = f 73 TT , < P r < 2tt , 

(4.1.2) 



e lP2 ° - 1 2tt „ 

U 2 S = _^ ,p o p" , < ^2 S < 27T . 

e 3 e*^ 2s — e 3 o 
Note that the Brillouin zone of the Q = 2 excitation is only the interval [^, 27r]. 

It is very difficult to evaluate the integrals in (4.1.1) explicitly except for a few special 
values of the momentum. One can check that 

P r = v r =0 =>- e r = 2|l-A], (4.1.3a) 

P 2s = v 2s =0 e 2s =4|l-A|, (4.1.3b) 

P r -> ^ =>- tv -> ± oo e r = 2(1 + A) , (4.1.3c) 

P 2s -> x v 2s ^oo e 2s = 4(1 + A) . (4.1.3d) 

This shows that even for infinite values of the rapidities the energy as a function of the 
physical variable momentum P is well-defined. Thus, as a function of P the Q = 1 
excitation has no singularity at P = y and also e 2s is continuous at both boundaries 
including P = For other values of the momentum than in (4.1.3) one can perform 
numerical integration where special attention has to be paid to the singularities in the 
integrand at t = 1. 

Proposition I: ([41]) All excitations above the groundstate can be composed in the large 
chain limit of the fundamental excitations Ei(P) := e r (P) and E 2 (P) := e 2s (.P) according 
to the following rules: 

m r m r m r 

AE Q)r (P, <j),ip) = J2 Eqw (P {k) ) , P = J2 pik) mod2n i Q = J2 Q(k) mod n - 
k=l k=l fc=l 

(4-1-4) 

subject to the Pauli principle, i.e. Q^ 1 ' = implies P^ ^ P^\ ■ 

It was argued in [103] by examining the zeroes of the transfer matrix that the same state- 
ment should also apply to the integrable three states chain although here the relation 
between rapidities and momentum is not known. 

For a detailed comparison with numerical methods and the perturbative approach to be 
presented below see [69]. 



4.2. High-temperature expansions 

In order to generalize the quasiparticle picture of the previous section we first study the 
low lying levels in the spectrum of the chiral Potts model perturbatively in this section. In 
particular, for n = 3 and n = 4 we calculate the dispersion relations of the lowest excita- 
tions in the charge sectors Q ^ 0. Some first results in this direction have been presented 
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in [68] for the self-dual version of these models. In this section we derive higher orders and 
admit general <fi ^ (p. We also present some explicit results on the next excitations. 

In [97] high-temperature perturbation series were computed for the disorder operator (or 
magnetization m) and the first energy gap in the momentum zero sector of the superin- 
tegrable Z3-chiral Potts model leading to exact conjectures for both of them. After the 
superintegrable chiral Potts model had been generalized to general n [70] perturbation 
series for the ground state energy, energy gap in the momentum zero sector, magnetization 
and susceptibility of this superintegrable Z n -chiral Potts model were presented in [83] . At 
the same time elaborate expansions of the ground state energy and some excitations of the 
superintegrable chiral Potts model for n G {3,4,5} and in particular perturbation series 
for the order parameters with general n were calculated in [7]. First perturbative results 
for the energy gaps at more general values of the angles 0, <p were obtained in [68] where 
second order high-temperature expansions for the translationally invariant energy gaps in 
each charge sector of the general self-dual Z3- and Z4-chiral Potts models as well as a first 
order expansion for the dispersion relations for general n was presented. 

For arbitrary n, N the groundstate of the Hamiltonian (3.1.1) in the limit A — > is 
given by: 

|GS) :=|0...0) (4.2.1) 

provided that — ^ < <p < \. For n = 3 (4.2.1) will be the groundstate for the larger range 
—n < if < n and for n = 4 (4.2.1) is the groundstate for — ^ < (p < 

The first excited states at A = for Q > and arbitrary P are the states 

|| a «»p:=||Q0...0» P (4.2.2) 

in the range — § < 9? < §. According to the definition of the Hilbert space 7i in [91], the 
states (4.2.2) give rise to proper eigenstates in the limit of AH^- n \ Thus, the corresponding 
gaps AEq : q(P) belong to the point spectrum of AH^- n \ 

More generally, we wish to argue later on that the complete spectrum can be explained 
in terms of quasiparticles. At A = 0, a single-particle excitation corresponds to just one 
flipped spin (4.2.2). Due to the absence of interactions /c-particle states have k flipped spins 
at A = 0. Therefore, one can think in terms of states keeping in mind that for A > one has 
to take the interactions into account using perturbation theory. Note, however, that except 
for the single-particle states and certain two-particle states we cannot directly compute the 
perturbation series for multi-particle states and check that they have the correct properties. 
Nevertheless, such a picture may be suggestive and is in fact traditional (see e.g. [35] for 
a transfer-matrix approach to the spectrum of the Ising model in d dimensions). 

First, one can derive a first order expansion for the dispersion relations of the lowest 
excitations in the charge sectors < Q < n using (4.2.1) and (4.2.2) [68]: 

AE Q , (P, a k , <(>) = «* (1 " " Q ") ) " T-pflJ cos ( P " ( X " v) + ° {X2) - 

(4.2.3) 



20 



From (4.2.3) we can immediately read off that in the limit A 
dispersion relations -P m i n ,Q is located at 



the minimum of the 



min,Q 



1 - 



2Q 



n 



(4.2.4) 



which demonstrates that for (f> ^ parity is not conserved. 

Let us now specialize to the Z 3 -chain. We will use the abbreviations: 



:= cos (|) , 



C := cos 



n 



K:=l- 4C Z 



U := 711 - 4 



cos 



rcj) 
~3~ 



(4.2.5) 



For n = 3 we can calculate the groundstate energy per site eo which is defined by 
perturbatively: 
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A 7 + C(A 8 ). 



(4.2.6) 

For (p = (f> = ^ (4.2.6) reproduces the result of [83]. The orders of the free energy per 
site eo are independent of iV if iV > k. Comments on the accuracy of a truncation to 5th 
order of (4.2.6) can be found in [82,91]. 

,2 

In principle, one can derive a critical exponent a for the specific heat from a per- 
turbation expansion of the ground state energy eo- In fact, a has been estimated using a 
13th order expansion of eo in [97] for the self-dual case 4> = tp and for the superintegrable 
case (f> = <f = ^ the Ising-like form of the eigenvalues has been exploited to calculate even 
higher orders of eo in [7] . The results in [97, 7] indicate a = g independent of the angles 

We omit the results presented for the lowest translationally invariant gaps of the three 
states chain presented in [91]. For n = 3 and general P we obtain from the states (4.2.2) 
the perturbation expansion (4.2.7) below for the dispersion relation of the lowest Q = 1 
excitation and the lowest Q = 2 excitation is given by (4.2) 3 ). In [69,91] we restricted 



3 ) The fourth order was independently obtained by G. von Gehlen using the slightly 
more compact notation of [69]. 
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to a 3rd order version of (4.2.7) because of the complicatedness of the 4th order which we 
present here just as an illustration. It is also not difficult to compute the 5th order on a 
computer but we do not want to present this to the reader here either. 
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(4.2.7) 

Note that the agreement even of the truncations to 3rd order of (4.2.7) with the results of 
a numerical diagonalization is usually good as was discussed in detail in [69] and in [91] 
for = ip = P = 0. 

We would like to mention that after specializing (4.2.7) to = ip = it becomes much 
simpler and it even makes sense to also present the 5th order [91]. For = <p = one 
can read off from (4.2.7) that for A small, a Klein-Gordon dispersion relation is a good 
approximation. It was shown in [2] that for = <p = and in the limit A = 1 one obtains 
£i(P) = 6| sin (|-)| which is also of the Klein-Gordon form on the lattice. However, one 
can show that for general A and = ip = a Klein-Gordon relation is nothing but a very 
good approximation to Ei(P) [91]. Also certain generalizations of Klein-Gordon dispersion 
relations can be ruled out [91]. 

In section 3.3 we mentioned that the kth orders of (4.2.5) - (4.2.7) are independent of N 
if iV > k + 2. In particular, this implies the existence of the limits iV — > oo of (4.2.6) and 
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(4.2.7) if the perturbation series converge at all. 

In the derivation of (4.2.7) we have not assumed that the Hamiltonian (3.1.1) is hermitean. 
Thus, we may admit e C. We have checked in a few cases that results of a numerical 
diagonalization at N = 12 sites are still in good agreement with (4.2.7) also for complex 
0. 

For the Z 4 -chain one obtains up to second order in A: 
AE h0 (P, 0, y>) = 2 (l + 2 sin {^^j ) " A2v/2 cos (p - 

2y/2cos (p+f)+3 CQS {2P -<!>)- 2 { , ,. 

_<1 4cos(^) + i + v ^ cos (f) lv'2cos(^) ' 

A£ 2 ,o(P,0, ¥>) =4v / 2cos(|) — A2 cos (P) 

_ 2 [ cos(P-0) + l cos(P + 0) + l _ 2 cos(2P)-l ] 3 

\l-v^sin(f) l + v^sin(f) l + V2cos(f) 4 v / 2cos(|)J 1 

(4.2.8) 

Again, AP 3;0 is given by (4.2). The special case <j> = (p and P = of (4.2.8) was presented 
in [68]. 

Returning to the Z3-chain we observe that for = <p = ^ we have to use degenerate 
perturbation theory. This was not done before [91] but can easily implemented for P = 
with the method of section 3.3 and was carried out in [91]. It is straightforward to verify 
the exact result (4.1.3b) which was first shown by [21] AP 2 ,o (f , f ) = 4(1 - A) up to order 
A 8 . Unfortunately, P 7^ is not yet accessible to perturbative computations because the 
lowest eigenvector of the potential has a very complicated P-dependence and becomes quite 
simple just for P = 0. Nevertheless, the fact that we are able to compute fairly high orders 
even in degenerate cases demonstrates the universality of the approach to perturbation 
expansions outlined in section 3.3. 

Also for the higher excitations we must apply degenerate perturbation theory. The 
next simplest case are the states where two spins are different from zero. For general n, 
P, — I < (f < ^ the space of the excitation with one spin flipped into charge sector Q 1 
and another one flipped into charge sector Q 2 is spanned by the states 

l|f? 1,ga »p:=IIQi ^0 Q 2 0...0» P , l<j<{"^' ifXtt (4 ' 2 - 9) 

j — 1 times 

Obviously, we will have to consider two cases: Qi 7^ Q2 and Qi = Qi- 

For Qi 7^ Q2 it is not obvious even for n = 3 how to diagonalize V in the space (4.2.9) in 
closed form although one can easily diagonalize it numerically for comparably long chains 
- for details see appendix A. 

In the second case Q\ = Q2 one can diagonalize the potential V in the space (4.2.9) 
explicitly exploiting a connection to graph theory (see e.g. [80]). We shift the details of 
the computation to appendix A and now present just the result. 
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It turns out that in the case of two identical excitations one has to distinguish between 
even and odd momenta in terms of lattice sites. It is convenient to introduce a further 
abbreviation dp encoding this distinction: 

PN PN 
dp 1 : = , if odd; 5% := 1 , if even. (4.2.10) 

27T 27T 

The result of the calculation in appendix A for the eigenvectors of the potential V is: 

2 / ,,„/ / ( 2 *- ^ ?)J' ^ ^„-i*o-l),, ^ Q, 



p 



H = sin 

V3 +(-!)" 



(4.2.11) 

The final result for the first order expansion of the energy for these excitations is for N > 3: 



(4.2.12) 



/ ^ (=?) 

+ !<*<[^±fzl" 

For a few remarks on the second order see the end of appendix A. 

4.3. Spectrum in the low-temperature phase 

In this section we calculate the ground state energy and the lowest excitations in the low- 
temperature phase using perturbative expansions according to [23] around A = of the 
Hamiltonian (3.1.2). We restrict once again to the Z3-version of the chiral Potts model 
(3.1.2) with periodic boundary conditions R = but impose no restrictions on the angles 0, 
cp. We present low-temperature expansions for the ground state energy and in particular 
the first translationally invariant energy gaps that before [82] had not been treated by 
perturbative methods because of high degeneracies. 

In each charge sector Q of the low-temperature phase there is one unique ground state. 
For arbitrary n it is given by (3.2.2) provided that — ^ < <fi < f • F° r n = 3 (3.2.2) is the 
ground state if — n < (p < n and for n = 4 (3.2.2) is the ground state for — ^ < </> < 
The excited states are more complicated and highly degenerate. The space of the first 
excitation is spanned by those states which have precisely two blocks of different spins. 
Furthermore, the values of the spins in these two blocks must have difference one. For 
fixed P, Q and — § < ^ < f we can choose the following basis for the space of the first 
excitation: 



\ a k 



1 n— 1 

Q ) := -= y^u l - Q \\ (l + l mod n)...(l + l mod n)l...l)) P . (4.3.1) 
V^^ ^ ' 



k times 
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In order to perform explicit calculations we now specialize to n = 3 with P = 0. In [82] a 
5th order low-temperature expansion of the groundstate energy per site eo was computed. 
Upon interchanging <p and ip it agrees with (4.2.6) as predicted by duality (theorem V). 
It is noteworthy that the expansion in powers of A of eo does not depend on the charge 
sector for large N. More precisely, the order A fc e fc ' ) of the free energy does not depend on 
the charge if N > k. However, ground state level crossings have been observed in [105, 66] 
which in the perturbative approach are due to the fact that for short chain length N high 
orders X k (k > N) do depend on the charge Q. These level crossings at fixed A in the 
massive phase were exploited in [66] to compute a critical exponent of the wave vectors by 
the argument introduced in [67]. For more details on the relation to perturbation series 
compare [82]. 

The calculation of the smallest gap AEq^ is more difficult. Let q be the projector 
onto the space spanned by the states | a^) (4.3.1). In this space, the potential acts as 
follows: 

qr(V)\a?) = -^(2e^ \cg) + e*J* \a%_ x )) 

gf(y) \a2) =-A( 2e ¥ | a? +1 > + 2e"^ (a^)) 1< k < N - 1 (4.3.2) 

Qf(V) lag.!) = -^=(2e"* \a%_ 2 ) + e~^ |a?». 
In the limit N — > oo the eigenvector with lowest eigenvalue converges to 

N-l 



k=l 

Using (4.3.2) and (4.3.3) we can calculate that for iV — > oo and < f 
w lta A Be , 1= 4V3co S (|)-A^cos(|) 

8cos(f)(l+co S (f)) 4(2 _ cos(¥)) } (1-3-4) 
3V3(3-4co S (|) ) 3V3co S (f) 



Comparing (4.3.4) with the corresponding high-temperature expansion (4.2.7) shows that 
up to the order calculated it coincides with A£ , 1)0 (A) + A£ , 2) o(A) at the dual point A = A 
with <fi, ip interchanged. Furthermore, one can argue that the finite-size corrections to 
(4.3.4) are of order iV -2 [82]. This further supports the identification of (4.3.4) with the 
dual of a two-particle state. Of course, (4.3.4) holds only for < ^. 

For > the states (4.3.1) are not the first excited states any more. Now we have to 
consider the following states: 



^( ||l...l 2...2...0...0)) P + w Q ||2...2...1...1)) 
+ . . . + c^" 1 ) 1 1 . . . . . . (n - 1) . . . (n - 1))) p^j . 



p 



(4.3.5) 
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Going through the same steps as before we find: 



lim AE Qtl =12 sin ( ^— ^ ) - W3cos 

N^co ' \ O J \ 6 / 



2 2 cos(f)-2 cos(f) + l l , 



7T 



(4.3.6) 



- A — < v + V- ,n ) + °( x )' 2 " < * 



V5\ cos (I) cos(V) 



(4.3.6) coincides with 3AEi j0 (\) at the dual point A = A with interchanged 0, </? in the 
corresponding high-temperature expansion (4.2.7) up to the order calculated. 

For = f the states (4.3.1) and (4.3.5) are degenerate. However, for large N the dominant 
contribution comes from the states (4.3.5) such that (4.3.6) is valid for = ^ as well. At 
(f) = (p = ^ this is in agreement with the exact result of [21]: 

lim AEq i = 6(1 - A) for <j> = p = J. (4.3.7) 



4.4. Evidence for quasiparticle spectrum 



Due to the studies of [41] (see proposition I) and [103] one expects a quasiparticle picture 
for the integrable Z 3 -chiral Potts model. Extensive numerical studies [68, 66, 69] showed 
that this quasiparticle picture is neither linked to integrability nor restricted to the Z 3 - 
chain. Furthermore, we were able to argue in [91] that this quasiparticle picture is to be 
expected in a quite general situation. Here, we do not want to go into technical details. 
We first recall the main statement, then discuss its origin and its consequences. 

By a Z n -quasiparticle spectrum we understand the obvious generalization of (4.1.4): 

Def. V: (Z n -quasiparticle spectrum) A family of spin chain Hamiltonians has a Z n - 
quasiparticle spectrum iff in the limit iV — > oo the spectrum of the excitation operator (3.1.14) 
can be obtained by the following rules: 

m r m r m r 

AE Q , r (P, ( j>, l p) = Y2E k (pW), P = J2P {h) mod27r, Q = mod n (4.4.1) 

k=l fe=l k=l 

with m fundamental quasiparticles carrying energy Ei(P) to E m (P) and charge to Q^ m \ 

Now the precise statement is: 

Proposition II: Let Hn be a nearest neighbour interaction spin chain Hamiltonian of the 
form Hn = H + XV such that Hn conserves Z n - charge for all A. Assume that H has a 
Z n - quasiparticle spectrum with n — 1 fundamental particles of charge = j (j ^ 0) and 
Brillouin zones [0, 2tt] . Assume furthermore that the perturbation series of the fundamental 
particles converge for A < Ao . 

Then for A < Ao the limit N — > oo of the complete excitation spectrum of Hn exists and it 
is a Z n - quasiparticle spectrum with n — 1 fundamental particles of charge = j (j ^ 0) 
and Brillouin zones [0, 2tt] . ■ 
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Remarks: 

1) For the Z n -chiral Potts model (3.1.1) with the parameterization (3.1.12) the assump- 
tions of proposition II can be easily checked except for a non-trivial radius of conver- 
gence which we shall check for the Z^-chain below in section 4.5. However, proposition 
II also covers more general cases, e.g. also the Hamiltonians (2.4.5), (2.4.6) and (2.4.9). 

2) Due to the assumptions the only overlap between proposition I and proposition II is 
the trivial point A = on the superintegrable line of the Z 3 -chain. 

Recall that one may expect to interpret the energy bands as continuous spectrum in the 
weak limit of the Hamiltonian [69, 91] and that the single-particle excitations (4.2.2) lead 
to point spectrum. The normalization factors for the two-particle states in (4.2.11) 
demonstrate that these states tend to zero for N — > oo and will therefore not give rise 
to proper eigenvectors. Furthermore, comparing (4.2.12) with the first order expansion 
for the single-particle states (4.2.3) one observes that this first order expansion of the 
two-particle excitations is in agreement with the quasiparticle rule (4.4.1). Up to first 
order the composite particle states satisfy either 2AEq$(P, <fi, (p) < AE 2 q^{2P, <p, ip) < 
2AE Q:0 (P + 2n, (f>,ip) or 2AE Qfi (P, </>, <p) > AE 2Q , k (2P, 0, <p) > 2AE Q:0 (P + 2tt, <j>, <p) 
depending on which one of the single-particle energies is larger. Thus, the two-particle 
states do indeed lie inside the energy band of two single-particle states and the boundaries 
are not included. Even more, we can see from (4.2.12) that the two-particle states become 
dense in this energy band for N — > oo. 

The proof of proposition II relies on the interaction in the Hamiltonian (3.1.1) be- 
ing very short ranged - in fact, only among nearest neighbours. In the massive high- 
temperature phase there is no spontaneous order and the correlation length is finite. Thus, 
if one puts two excitations of 'short' chains with a sufficient separation on a longer chain, 
the interaction will be negligible. In pictures: 

N M N + M 



- A , _ , * /v - 



•- 



EuPu Qi E 2 , P 2 , Q 2 n,m^co El + E 2 , P x + P 2 , Q 1 + Q 2 

(4.4.2) 

For example, putting one single-particle excitation one the left half of the chain and another 
on the right half will approximate a two-particle excitation. 

In order to make the derivation rigorous one has to show the vanishing of boundary terms 
in (4.4.2). This can indeed be done using e.g. perturbative arguments [91]. Because this 
is very technical we do not want to recall the complete argument here. Just note that the 
crucial point is that the momentum eigenstates have normalization factors iV", M~%. 
Any operator acting only at boundaries yields only a finite part of these states in contrast 
to the operators T N and AH^' which act on the complete chain and yield complete 
momentum eigenstates. The finite pieces of momentum eigenstates are suppressed by the 
normalization factors N~? in the infinite chain length limit. 

Note that boundary terms are substantial for conformally invariant systems with long 
ranged correlations. Also in the low-temperature phase boundary terms play an important 
role because the free part of the Hamiltonian depends on the difference of neighbouring 
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spins [82]. Thus, proposition II applies neither to critical points where one might have 
conformal invariance nor to the low-temperature phase. 

We have not assumed the Hamiltonian to be hermitean. In particular, the quasiparticle 
picture should also be valid for e C as long as the single-particle excitations exist and 
converge. This is indeed supported by numerical calculations [126]. 

The argument proving the quasiparticle structure can be refined in order to give an 
upper estimate for the rate of convergence in N of the energy of a /c-particle state. As an 
approximation to a /c-particle state for kN sites, total energy E tot and total momentum P 
we may take the fc-fold tensor product of single-particle states 

||fciV; £ t ot» P := \\N; E 1 }} Pl <g> . . . <g> \\N; E k )) Pk (4.4.3) 

with E tot = Xw=i Eh P = Xw=i Pi- Now, the deviation from the limit N — > oo is given 
by: 

k 

P ((kN;E tot \\ AH$\\kN;E tot )) P - E tot = J] Pl ((AT; E t \\ G(AH N )\\N; E t )) Pl = 0{N~ k ) 

i=i 

(4.4.4) 

where O(AHn) is some operator that acts only at sites 1 and N. The first equality simply 
uses the definition of the scalar product in tensor products. The last equality is more 
profound and due to the fact that operators acting only at boundaries of the chain are 
suppressed by iV -1 due to the normalization factor in the finite fourier transformation for 
momentum eigenstates. This general argument is confirmed by our results for the two- 
particle states. Expanding cos(x) = 1 — ^x 2 + 0(x 4 ) we can read off from (4.2.12) that 
the first order correction of the rth two-particle state with respect to the boundary of the 
energy band behaves as N~ 2 which agrees with (4.4.4). 

In summary, we have the following results on the finite-size effects: 

1) Single-particle states converge exponentially in N as long as the corresponding per- 
turbation series converge since the rth orders are independent of N if N > r + 1. 

2) The deviation of the energy of a /c-particle state (k > 1) from the limit is at most of 
order N~ k for N — > oo. Note that consequently, the iV-dependence of some energy 
eigenvalue yields only a lower bound on the number k of particles involved. 

In particular, the energy of any state remains unchanged to order ^ in a finite-size 
system. Thus, the only modification in (4.4.1) at order in the massive high tempera- 
ture phase is a discretization of the momentum (and possible minor modifications of the 
Brillouin zones and selection rules [103]). 

Note that the proof of the vanishing of boundary terms as sketched above and pre- 
sented in detail in [91] also directly applies to the Hamiltonian (3.1.1) itself. So far, we 
have restricted ourselves to periodic boundary conditions but one could also use cyclic, 
twisted or even free boundary conditions. Our argument shows that all these different 
choices lead to the same spectrum in the limit iV — > oo. In particular, our results are 
valid for all choices of boundary conditions and one is free to choose those which seem 
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most appropriate, e.g. one can leave the ends of the chain open instead of the unnatural 
end-identification for a realistic physical system. 

Again, this observation for the massive high-temperature phase is to be contrasted with 
other situations. In particular, at the second order phase transition (p = <p = 0, A = 1 the 
correlation length becomes infinite and the boundary terms are very important [37, 74, 90]. 
Even in the massive low-temperature phase one observes long range order and boundary 
terms cannot be neglected [82]. 

Similarly one can directly argue that in the low-temperature regime all charge sectors are 
degenerate for N — > oo, at least if the perturbation expansion converges [82]. That means 
that in the massive low-temperature regime the thermodynamic limit does not depend on 
Q. The independence of the spectrum of the charge Q in the low-temperature regime is 
also expected due to duality (theorem V) because in the massive high-temperature phase 
it does not depend on the boundary conditions R. 

So far we have not addressed the question of whether the fundamental particles satisfy 
a Pauli principle or not - note that the above discussion is intrinsically insensitive to a 
Pauli principle because the limit was defined such that the spectrum forms a closed set. 
Nevertheless, for the special case n = 3 and = ip = ^, eq. (4.4.1) was obtained in [41] 
supplemented with the Pauli principle. Fortunately, due to (4.2.12), we have some control 
over the finite-size dependence of the scattering states of two identical particles in the 
general case. Up to first order in A these finite-size effects do essentially neither depend on 
the charge Q nor on the number of states n for periodic boundary conditions. Therefore, 
the nature of the fundamental excitations can be determined by looking at one particular 
choice of Q and n. However, for n = 2 one obtains the Ising model where it is well-known 
that the excitation spectrum can be explained in terms of one fundamental fermion (see 
e.g. [104]). This indicates that the fundamental excitations for general n should behave 
similar to fermions. In particular, for a scattering state of two identical excitations i and 
j the momenta must satisfy Pi ^ Pj. In a scattering state of two different fundamental 
particles these two fundamental particles can easily be distinguished because they carry 
different Z n -charges. Therefore, two different particles should not be subject to a Pauli 
principle (like it is the case for two different non-interacting fermions). 

From theorem V we conclude that the quasiparticle interpretation for the high-tem- 
perature phase of the general Z n -chiral Potts quantum chain given by proposition II can be 
pulled over to the low-temperature phase. The duality transformation interchanges charge 
sector Q and boundary conditions R. Thus, the ground state of the high-temperature phase 
is mapped to periodic boundary conditions R = in the low-temperature phase. However, 
the fundamental excitations are mapped to different boundary conditions corresponding 
to the charge sectors R G {1, . . .n — 1}. Therefore we observed only composite particle 
states in section 4.3. 

4.5. Convergence of single-particle excitations 

As far as the quasiparticle picture of proposition II is concerned the main open question 
is the convergence of the single-particle states, or equivalently the existence of the limits 
N — > oo of the corresponding eigenvalues of the Hamiltonian. We have argued in section 4.2 



29 



that convergence of the perturbation expansions is sufficient to guarantee the existence of 
the limits N — > oo. Therefore we will discuss the radius of convergence for the perturbation 
expansion of the single-particle excitations in this section. 

The potential for AH^ as defined in (3.1.1) and (3.1.14) is unbounded if N is not fixed. 
Thus, we have to apply the Kato-Rellich theory of regular perturbations. Reviews of this 
subject can be found e.g. in the monographs [117,100]. The main results we are going 
to use were originally published in [118,99]. The theory of Kato and Rellich applies in 
particular to operators of the form H(X) = H + XV. 

Suppose that the single-particle eigenvalues AE have a non-zero distance from the 
scattering eigenvalues (the continuous spectrum) at A = 0. Then it is clear from the 
discussion in the previous sections that these eigenvalues are non-degenerate and isolated. 
In particular, the resolvent (AH^\X) — z)~ x is bounded for \AE — z\ > 0. Restricting 

to the hermitean case, this is sufficient to guarantee that the AH^\\) are an analytic 
family in the sense of Kato. In this case, the Kato-Rellich theorem ([117] Theorem XII. 8) 
may be used to guarantee a non-zero radius of convergence ro > for the single-particle 
eigenvalues of AHj^'(X). 

In order to obtain explicit estimates of the radius of convergence one needs the inequality 

|| V | a) || < V\\H | a) || +/C|| | a) || (4.5.1) 

on T>(Hq) which in our case is dense in the the complete Hilbert space 7i. Then, the 
isolated point eigenvalues of H(X) are convergent at least for 

A<n:=V- 1 (4.5.2) 

as long as these eigenvalues do not come in contact with continuous spectrum [99]. On 
the one hand this criterion is very simple, on the other hand one must estimate not only 
the constant V but also examine the level crossings between single-particle excitations and 
scattering states. There is another estimate r 2 that guarantees the separation of eigenvalues 
as well but gives smaller radii of convergence. For self-adjoint H$ with isolated eigenvalue 
£q°' ) where the nearest eigenvalue E± has distance e := \e[ ^ — Eq°^\ (e _1 = ^(E^) ||) 
the perturbation expansion of Eq(X) is convergent for 



A < r 2 : = 




(4.5.3) 



and there are no crossings with neighbouring levels. In order to compare the estimates 
(4.5.2) and (4.5.3) let us assume /C = and l-E^I = e - For this almost optimal case one 
has r\ = 4r2 showing that the criterion (4.5.3) is much more restrictive. 

Let us now apply these general results to the present case of Z n -spin quantum chains. 
For non-degenerate single-particle eigenvalues the Kato-Rellich theorem can be applied to 
guarantee a positive radius of convergence ro- Then we know from section 4.4 that the 
spectrum of AH^\X) is a quasiparticle spectrum for A < r . This fact can be used to 
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calculate the constant V and obtain explicit estimates r\ (where level crossings still have 
to be discussed) or r 2 . One can obtain the estimate (4.5.1) with /C = using Schwarz' 
inequality: 

V := sup <°' A V' a) ■ (4.5.4) 
W6« IIA^b I a ) II 

In general, this supremum need not be finite but then it is very difficult to ensure con- 
vergence at all. In our case, the important observation is that due to the quasiparticle 
picture we can evaluate (4.5.4) exclusively from the single-particle excitations. To see this 
one performs a first order expansion in A for any composite particle state, compares co- 
efficients and uses the quasiparticle property to expand the expectation values of AH^\ 
and AV in single-particle excitations. Thus, V can be calculated as 

V-^^ljWK (4.5.5) 

In order to implement this program explicitly we specialize to the case of Z 3 with 
the parameterization (3.1.12). At A = both single-particle eigenvalues are isolated for 
— 77 < if < ^. This guarantees a non-zero radius of convergence r . 

The simplest case is the parity conserving case <j> = ip = 0. Here, the maxima are located 
at zero momentum P = and both charge sectors are degenerate. Furthermore, we have 
\\AH$ \\s Q )) o|| = e = From (4.2.7) at <j> = cp = we can therefore read off V = §, 

or in terms of radii of convergence 

3 3 

fi = -, r 2 = -, for n = 3, = (p = 0. (4.5.6) 
Z o 

r 2 = 0.375 is certainly too small (for more details see [91]). 

For general angles < <p> < n, the free part of the Hamiltonian || AHjy\\\s®)) o|| is min- 
imized for Q = 1 and the potential p((s^|| AV||s^)) p is maximal for P = Thus, we 

read off from (4.2.7) V = (v^sin (^-^)) • The relevant distance e in the spectrum at 
A = is given by the distance from the Q — 2 fundamental particle to a state of k particles 
with Q = 1 where k depends on (p. Note that k must be equal 2 = Q mod 3. Thus, one 
has e = 4|A;sin - sin (£±£) | and E^ 0) = 4 sin (^). This amounts to the following 

radii 

r (*-<p\ v / 3sin(^)|A;sin(^) -sin(^)| 

r 1 = y/3sm[ , r 2 = min — , \ ' , 1 4 n 4 - 5 - 7 

V 3 7 ' fc 2 (sin(^) + |A;sin(^) -sin(^)|) 1 ' 

for n = 3, < <p < n. For < ip < ^ the closest state is a two-particle state, i.e. k = 2. 

For — > tj the situation is contrary to that at = ip = 0. The Q = 2 particle state 
becomes degenerate with two Q = 1 scattering states at <p = ^ such that the radius of 
convergence must tend to zero for <p — > 77. Whereas r 2 has precisely this property, r\ tends 
to 0.866 . . . which is certainly too large. 
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Because for small <p we would prefer the large radius of convergence r\ but at <p ~ ^ 
this is much too large and r 2 seems more appropriate we have to enhance the estimate 
given by r\ by a discussion of level crossings between single-particle states and scattering 
states. For < <p < ^ the first level crossing of this kind will take place between the Q = 2 
single-particle excitation and a two Q = 1 particles scattering state. 

It is very difficult to determine those values of A explicitly and precisely where they take 
place. Therefore, we will use the first order approximation of the perturbation expansion. 
We are looking for those values of A where a single point P exists such that x(k) := 
/cA£ , i j o(y, 4>i V 7 ) — A.E2fi(P,4>,ip) (with k = 2 mod 3) vanishes. The fact that we are 
looking for no real crossings but x(k) = implies = 0. Inserting (4.2.7) and (4.2) up 
to first order leads to the condition 



sin (¥-l)= sin ( P + l) (4 ' 5 ' 8) 

Eq. (4.5.8) has a solution 

P = P- (4.5.9) 

k+1 v ' 

that does not depend on 0. Now we can solve the linear equation x(k) \\ = for the value 
Ao- One obtains 

AoW = *m (4 . 5 . 10) 

™ s (ra + !)- fccos (ra-f) 

The k = 2 special case of (4.5.7) - (4.5.10) was discussed in [91]. 

Fig. 2 shows a plot of the estimates (4.5.7) and (4.5.10) for the self-dual case (p = (p. Note 
that r\ and r 2 are independent of (j). However, we have assumed that the Hamiltonian is 
hermitean and therefore must be real. 

For completeness we have also included an estimate for the boundary of the massive high- 
temperature phase in Fig. 2. At this boundary, levels of the Q = 1 particle with generically 
non-zero momentum cross with the ground state. Its explicit location has been obtained 
estimating the minimum of the dispersion relation (4.2.7) with P = ^ and solving the 

second order approximation A£'i j o(f , 4>, <j>) = for A. This is a slightly modified version 
of the estimate already proposed in [97] where we have used the second order because it 
gives good agreement with exactly known points. 

For a detailed discussion of Fig. 2 in the interval < ip < ^ see [91]. 

The level crossings transition Ao divides the massive high-temperature phase of the Z3- 
chiral Potts model into two parts which we label Ha and lib. Note that above the super- 
integrable line V 7 = f there are infinitely many 'windows' of type Ha in which the Q = 2 
quasiparticle exists in the entire Brillouin zone [0, 27r]. At those points where the Q = 2 
fundamental particle crosses k Q = 1 particle states (k = 2 mod 3) already for A = 0, the 
radii of convergence become zero. The first crossing of this kind takes places at ip = ^ 
with a two-particle state. The crossings and consequently also the 'windows' accumulate 
for (p — > n. Note that for (p > ^ the estimates r 2 and Ao follow each other closely, i.e. the 
estimate r 2 is already close to optimum. 
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Boundary of the phase 



Ila 



Ila 



r2 



lib 



Fi;£. 2: Radii of convergence and boundary^©! the ntasslve higrjj^ipperat'ur'e v 



ffih kse for the liemiitean A^-chain. — r\ is an estiMate ensuring convergeiice 

ncOlevel crossings betweeK/$oint and continuous ufiUctrum occur. The es3tiJ5i4te 
r 2 also ensures the absence of level crossings(f> =TJpe perturbation series are 
definitely convergent for A < r 2 although the true radius of convergence is 
larger. It extends until the value Ao where the first level crossings between 
fundamental quasiparticles and scattering states occur. 
Note that r\ and r 2 are independent of <fi whereas for Aq we put ip = 4>. 



ljce if 



'■/-■>'..\.V^,W.-> 



In part Ila the derivation of the quasiparticle picture is rigorous. Thus, in regime Ila 
the spectrum is a quasiparticle spectrum with two fundamental particles existing for all 
momenta. In [69] we have presented numerical evidence that regime lib probably also 
exhibits a quasiparticle spectrum with two fundamental particles where the Q = 2 particle 
has the unusual property that it exists only in a limited range of the momentum P. 
At (f = (j) = ^ this statement is given by proposition I. We expect that the idea to 
approximate multi-particle states by putting single-particle states of 'small' chains with 
a sufficient separation on a longer chain and to use the finite correlation length in order 
to ensure vanishing of boundary terms (which we cannot show directly like in [91] any 
more) will apply also in regime lib for general angles <p, (p. However, in contrast to section 
4.4 we loose control over the fundamental Q = 2 excitation because the perturbation 
series does not converge any more and there is no guarantee for the completeness of this 
construction. At least it is plausible to still expect a quasiparticle spectrum in regime lib 
with two fundamental particles of which the Q = 2 particle may have a Brillouin zone that 
is smaller than the interval [0, 27r]. 
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Tout cela est d'une grande beaute et d'une extreme importance, mais malheureusement nous ne le 
comprenons pas. Nous ne comprenons ni I'hypothese de Planck sur les vibrateurs, ni /'exclusion des 
orbits non stationnaires, et nous ne voyons pas, dans la theorie de Bohr, comment, en fin de compte, la 
lumiere est produite. Car, il faut bien I'avouer, la mecanique des quanta, la mecanique des discontinuites, 
doit encore etre faite. 

H.A. Lorentz in "L'ancienne et la nouvelle mecanique" (1925) 

5. Correlation functions 



In recent papers a systematic investigations of the correlation functions of the Z3-chiral 
Potts model in the massive phases has been started. First, a non- vanishing wave vector has 
been predicted in [66, 105] for the massive low-temperature phase and its critical exponent 
was calculated from level crossings. Below we wish to present the direct computations 
of the correlation functions of [82,69,91]. Note also that for the massless phases around 
A ~ 1 of the Z3-chain correlation functions have been derived in [3] using results from 
conformal field theory. 

5.1. Correlation functions of the high-temperature phase 

In this section we study correlation functions for the Z 3 -chiral Potts model perturbatively. 
Before defining correlation functions, we first note that the two-point functions are trans- 
lationally invariant because the groundstate \v) is translationally invariant: 

(v\ r+ +r r r \v) =(v \ r+ +1 r!|i;) 

(v\ a+ +r a r \v) = (v | \ v) . 

Thus, it makes sense to define the correlation function for an operator S by the following 
expression: 

= H5++i5iH _ W5^, |.) H5.M < x < — (5.1.2) 

{v\v} {v\v) z 2 

where | v) is the eigenvector of the Hamiltonian to lowest energy. Here, we do not assume 
that | v) is normalized to one and have therefore included the proper normalization factors 
in (5.1.2). 

The operator H for the Z 3 -chiral Potts model can be either r or a. (For n > 3 also different 
powers of these operators may be interesting.) The correlation functions of the operators 
T x and a x have the property 

Cr(-x)=C r (xy , 

c a (-x) = c a (xy = c a (x) 

such that it makes sense to restrict to positive x. Note that (5.1.3) follows by complex 
conjugation using (5.1.1). Explicit calculations show the validity of (5.1.1) and (5.1.3) as 
well. 
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(5.1.1) 



(5.1.3) 



Let us now turn to the explicit computation of correlation functions for the Z 3 -chain. 
In order to be able to calculate the correlation functions we need to know the groundstate 
| v) . We will calculate it from the free ground state | GS) using the perturbation expansion 
(3.3.5). We should stress again that although we assume the free groundstate | GS) to be 
normalized to 1, this is not necessarily true for the complete state \ v) . The expansion of 
the groundstate | v) provides us with an expansion for the correlation functions in powers 
of A 

oo 

C E (x) =J2^ci u) (x). (5.1.4) 

Note that according to (3.3.5) a kth order expansion of the groundstate yields a k + 1th 
order expansion of the groundstate energy as a byproduct. 

The operator T x creates charge such that charge conservation implies (v | T+ \v) = (v \ 
T x \v) =0 for all x. Thus 



M 


r+ 


■iri 


v) 


(v\ 


v) 



C r (x) = , , V • ( 5 - L5 ) 



Here we omit the intermediate steps and present just the final results of the compu- 
tations. The interested reader may find some more details in [91]. 

The final result for the expansion of C a (x) is 
Cf\x) =CW( X ) = 0, 

C?\ X ) =^ {4,o + <?} , C?\x) = ^ [s x , + 5 -f), 

r( 4), x 1 / c f 2(1 -IPC 2 ) 3\ { 1 + 20C 2 n (5.1.6) 

/ 2(1 + 2C 2 ) ; J_ 

+ X ,2 y + 16C2 

For the other correlation function Cr(x) one finds the following final result 

i± ( A 2 <t> 

Ml), > 1 f , it S 1 8 \ , _,± I 2 8 \ t be'* 



■lib , o x ■ 4<* • 2rb ■ 4ti 

8e l ^(19C 2 -4) 3e"^-20e 1 ^ 3e" 1 ^ 



+ 4,2 I — + 



m? 8C 2 n 

7^ ^ 4C~i J ^ Ux > 4 ~&p 



. ,'40C 2 -7 9 \ 35e^ 

+ X ,3 ^ + 77^2 + d M 



(5.1.7) 
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It is easy to see that C^\x) and C^\x) are independent oi N \i N > 2k and x < k. 

C a (x) is real and positive for all values of and y? up to the order calculated. However, 
it is not easy to read off from (5.1.6) what might be the form for large x. Thus, we specialize 
to 4> = 4> = 77 and calculate two further orders for C a (x): 

Cf\x) =C^\x) = Cf\x) = Cf\x) = , 
1 



c<?H 



X 



9 



{44,o + 4,i} 



1 



C^\x) =^-{54,0 + 24,2} , 



(5.1.8) 



1 



6561 



{1904,0 - 134,1 + 384,2 + 604,3} • 



Cr(x) in general has a non- vanishing imaginary part and therefore is worth while 
being considered in more detail. Thus, we specialize again to the superintegrable case 
4> = if = 5 4 ) and obtain after calculating two further orders 



4°\x 



4 2 \x 

C^(x 

Cfix 
c£\x 



C^(x) 



— 4,o , 



4 1 \x) = 5 x , 1 I | + z- 9 



1 .Vs^ 



{4,i + 25 x , 2 } + {-S Xtl + 24,2} , 

i {44,i + 105 x ,2} + {44,i - 104,2 + 204,3} , 

1 /3 
= — {274,i + 184,2 + 2104,3 - 704,4} + {-274,i + 184,2 + 704,4 , 



2187 



{454,i + 1084,2 + 2524,4 - 1264,5} 



+ i—— {154,i - 364,2 - 144,3 + 844,4 + 424,5} 



2187 



1 



39366 
. V3 



{3814,i + 2144,2 + 23144,3 + 7844,4 + 23104,5 - 18484,e} 



+ 1 



39366 



{-3814,i + 2144,2 - 7844,4 + 23104,s} • 



(5.1.9) 

Of course, we still have to calculate the sum (5.1.4). Thus, changes of signs in individual 
orders need not necessarily turn up in the final result. In fact, it turns out that the 
imaginary part of Cr(x) is always positive up to order 6 because the smallest orders are 
positive and they dominate the others. However, for sufficiently small A the real part does 
indeed change signs around x = 4. We fit (5.1.9) by a complex exponential function: 

C r (x) = a e^~^) x + (1 - a)4,o , (5.1.10a) 
C a (x) =p e'-t +qS xfi . (5.1.10b) 



4 ) Another natural specialization would be the standard Potts case (j) = = 0. This was 
discussed in detail in [91]. There, a deviation from the well-known relation (see e.g. [104]) 
£r ~ m(A) -1 was observed which is clearly due to lattice artifacts. 
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In (5.1.10) we have also taken into account that from (5.1.8) c " \i) ~ ^ independent of the 
correlation length 

If (5.1.10a) is the correct form for Cr(x) we infer from (5.1.9) that L is about 14 for small 
A. We can also see from the higher orders that L increases with increasing A such that it 
might well be singular at A = 1. The correlation length £ tends to zero as A — > 0. This 
implies that - after proper re-normalization of the Hamiltonian - the mass gap becomes 
infinite at A = 0. It has already been observed in [68] that there are physical reasons to 
divide (3.1.1) by v^A which would have exactly the effect of infinite mass at A = 0. Fits to 
(5.1.10) for A G {j, |, | } in the superintegrable case are given by the values in table 1. 



A 




a 


L 


to 


V 


Q 


1 mm 


L ^rn i n 

2tt 


0.25 


0.55(3) 


0.55(5) 


14.3(2) 


0.25(2) 


0.35(4) 


0.32(4) 


0.471 


1.07(2) 


0.50 


0.9(1) 


0.59(3) 


16.5(8) 


0.38(4) 


0.35(3) 


0.24(3) 


0.401 


1.05(5) 


0.75 


1.5(6) 


0.64(3) 


18.3(8) 


0.55(6) 


0.36(2) 


0.09(2) 


0.308 


0.90(4) 



Table 1: Parameters for the correlation functions (5.1.18) at = cp = \. 



The estimates in table 1 have been obtained as follows. First, £r has been estimated by 
calculating Re ( In ( qJ^+ x ) ) ) ~ 1 and averaging over x. Next, the zero of Re(e«r Cr(x)) has 
been estimated by linear interpolation for two neighbouring values and L/4 was obtained 
by averaging. Finally, a was estimated such that the difference 



Ee(C r (x)) - ae~t cos (^j^ (5.1.11) 

is minimal for x = 1,2. That this procedure yields reasonable fits is demonstrated by 
Fig. 3 which shows the stretched correlation function e^r Cr(x) in comparison to the fits. 

6 — x 

The 'error bars' are not really error bars but given by ae~^ r which gives an idea how much 
the values have actually been stretched and what might be the contribution of the next 
orders in the perturbation expansion. The agreement for all x not only in the real part 
but also in the imaginary part is convincing. 

Let us now discuss the implications of (5.1.7) under the assumption that (5.1.10a) is the 
correct form for general values of the chiral angles. From the leading orders in (5.1.7) we 
read off the following identity for the ratio of Ct(1) and Cp(2): 



Cr(2) = gA 2 + 0(A 3 ) = e4 x , , om2 
C r (l) e4\,n(\2\ 2C 



X + 0(X z ). (5.1.12a) 



3C 



X + 0(X 2 ) 



On the other hand we immediately obtain from (5.1.10a) 



Cr(2) 
Cr(l) 



. J_ gjri 

€r e L 
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(5.1.12b) 
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exp(x/^) Re(C r (x)) HO— I 

exp(x/^) Im(C r (x)) H — I 

0.59 cos(2 71 x/16.4) 

0.59 sin(2 71 x/16.4) ----- -■ ^ 

Fig. 3: Correlation function Cruel stretched by in comparison to the 
J 1 ^ 1 L 



m (5.1.10a) at ^ = ip = f , A = f. ^ 
^which convey! an idea how^nuch the va lues have actuary been stretched. The 



he "error bars' are given by ae «r 
)S have a< 

oscillatory contribution to Cr(x) is clearly visibfe. 



<> 



Comparison of (5.1.12a) and (5.1.12b) leads to 

6tt 



L = 



Re(0) 



In 



2cos(f) 



lm(0) 
3 



(5.1.13) 



for small values of A. It is noteworthy that we obtain the same result for the oscillation 
length if we apply a similar argument to i n lowest non- vanishing order with x\,x<l G 

{1,2,3,4}. At = f (5.1.13) yields the approximations L = 12, £ r = 0.52, 0.80, 1.2 for 
A = 0.25, 0.50, 0.75. The agreement with the numbers of table 1 is very good. Thus, 
for very high temperatures the oscillation length L is proportional to the inverse chiral 
angle </> _1 . In particular, the oscillation vanishes smoothly for — > 0. According to (4.2.4) 
for very high temperatures the minimum of the dispersion relation of the fundamental 
particles is also proportional to (p - for n = 3 and Q = 1 we have P m i n = Re ^ . Thus, we 
obtain from (5.1.13) for very high temperatures 



A^0 r 



2tt 



V0, (p. 



(5.1.14) 



Furthermore, the second order in (4.2.7) shows that the minimal momentum P min decreases 
with increasing A (compare also [69]). Similarly, we read off from (5.1.7) that the inverse 
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oscillation length L -1 also decreases with increasing inverse temperature A. Thus, (5.1.14) 
has a chance to be valid for all values of A in the massive high-temperature phase. Indeed, 
using the values of -P m i n given in table 8 of [68] we see that P m \ n L = 2w holds quite 
accurately for A = 0.25,0.5,0.75 at = <p = ^ (compare table 1). We suspect that in 
general this relation is not exact but an excellent approximation. 

Note that even at = ip = ^ the correlation lengths £r and £ CT are clearly different. 
Furthermore, £ CT coincides with its dual in the low-temperature phase whereas £r does 
not (see [82]). This can be understood in terms of the form factor decomposition to be 
presented in section 5.4 (see [91]). 

5.2. Numerical computation of the correlation functions 

In this section we study the correlation function for the operator V of the Z 3 -chain numer- 
ically. We will in particular check the relation (5.1.14) numerically. 

In order to evaluate the correlation function (5.1.5) numerically one needs the groundstate 
of the Hamiltonian (3.1.1). The groundstate of the Z3-chain is easily obtained for general 
values of the parameters (even in the non-hermitian case cj) e C) using vector iteration up 
to N = 13 sites. The more difficult point is to calculate the matrix elements of r^^Ti 
because this operator does not conserve momentum and thus does not leave a space of 
momentum eigenstates invariant. It is crucial to really work with r^ +i ri and not to 
replace it by X)r=i ^t+r^r because this would destroy the oscillation. Table 2 shows 
the correlation function Cr(x) obtained in this manner for N = 12 and N = 13 sites at 
</> = </?= : |,A=i which corresponds to Fig. 3. This table also contains the perturbative 
results for the correlation function Cp 6rt ' (x) which were obtained in section 5.1. 



X 





1 


2 


3 


4 


5 


6 


C* ert (x) 


1 


.18868+. 07384i 


.04561 + . 03980i 


.00992+. 01747i 


.00091 + . 00674i 


-.00088+. 00263i 


-.00074 


Cf' num -(x) 


1 


.18881 + . 07385z 


.04587+. 03967i 


.01004+. 01737i 


.00126+. 00679i 


-.00056+. 00224i 


-.00080 


C^ 3 ' num '(x) 


1 


.18882+. 07385i 


.04588+. 03967i 


.01007+. 01738i 


.00132+. 00684i 


-.00043+. 00242i 


-.00063+. 00058i 



Table 2: Perturbative results and numerical results for iV = 12 and N = 13 sites for the correlation 
function Cr(x) at = cp = ^, A = 



The agreement between the results of both methods and also between iV = 12 and iV = 13 
sites is good. This shows that on the one hand higher orders are indeed negligible in (5.1.9) 
for x < 7 and on the other hand that the finite chain length does not considerably affect 
the correlation function Cr(x). Thus, one can take the values obtained e.g. at iV = 12 
sites as an approximation for the infinite chain limit as long as x < 6. 

Table 3 collects results of numerical calculations at iV = 12 sites for the self-dual values 
<j> = (p E {x> !> T"} an< ^ a ^ ew integrable cases. 
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P . 

1 mm 


L ^min 

2tt 


3tt 
4 


3tt 
4 


1 

4 


0.62(8) 


10(1) 


0.4(2) 


0.746 


1.18 ±0.12 


3tt 
4 


3tt 
4 


1 

2 


1.2(6) 


12(1) 


0.4(1) 


0.685t 


1.34 ±0.14 


3tt 
4 


3tt 
4 


3 
4 


1.4(1) 


16(1) 


0.65(3) 


0.602t 


1.54(8) 


7T 

2 

7T 

2 

7T 

2 


7T 

2 

7T 

2 

7T 

2 


1 

4 
1 
2 
3 
4 


0.6(1) 
0.60(5) 
1.3(2) 


14.05(6) 
16(1) 

29(6)* 


0.44(8) 
0.59(3) 


0.470 
0.401 
0.189 


1.051(4) 
1.07(8) 
0.9(2) 


3tt 
8 

Q _ 

8 

3-7T 

8 


3?r 
8 

o7T 

8 

3tt 

8 


1 

4 

1 
1 

2 
3 
4 


0.53(3) 
0.8(2) 
1.8(6) 


18(2) 

25(1)* 

c;n -i- oq* 
oU ± zo 


0.62(1) 


0.341 

0.283 
noil 


1.0(1) 
1.14(6) 
1.7(8) 


3tt 
8 


-% 0.98942 


1 

4 


0.7(2) 


oo 









3tt 
8 


0.69919 


1 

2 


1.1(3) 


49(7)* 




0.159 


1.2(2) 


3tt 
4 


7T + z 1.70004 


1 

4 


1.0(5) 


5.9(1) 


0.3(1) 


7T 

3 


0.99(2) 



Table 3: Parameters for the correlation function (5.1.10a) 
calculated numerically for N = 12 sites. 



The real part of the correlation function Cr has no zero in the interval [0,6] for the entries 
marked with a '*'. Therefore, the procedure of section 5.1 could not be used for the 
computation of L. Here L has instead been obtained by the formula 

L = -Y — ^. (5.2.1) 

The values of -P m i n in table 3 have been obtained by first calculating Ei(P) numerically 
at N = 12 sites and afterwards minimizing the finite fourier decomposition of (Ei(P)) 
numerically. This method for estimating P m i n failed for those values marked with a 't' 
because here Ei(P) partly crosses two-particle scattering states. Ei(P) becomes even 
negative for = cp = A = |. This implies that the groundstate does not belong 
to the Q = P = sector and this point in the parameter space lies in the massless 
incommensurate phase. However, we have ignored this fact and used the lowest energy 
state in the Q = P = sector for evaluation of the correlation functions. This is justified 
because it yields results that are similar to those for the other values of the parameters. The 
values of P m ; n marked with a 't' have been determined by first determining the smallest 
energy gap in the Q = 1 sector at finite N = 8, ... , 12. For 8 < iV < 12 the minimum of 
the discretized dispersion relation is located at P = Next, a polynomial interpolation 
between these five values for the energy gap was minimized numerically. 

That the above procedures yield reasonable fits was demonstrated by a figure in [69] . 

We have also included a few values for integrable points at the end of table 3, i.e. points 
where the parameters satisfy cosy? = Acos</>. For ip = ■^ L , A = -|, the integrable case is 
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hermitean and behaves precisely like the self-dual case <f> = (p. For small values of A one 
has to choose one of the angles ip, <fi complex in order to satisfy the integrability-condition 
(3.1.13). We choose <p complex: For p < ^ we choose purely complex <p £ iR whereas for 
<p > \ we take ifr — irEiM (compare also [126]). For complex choices of the Hamiltonian 
becomes non-hermitean such that the excitation energies become complex. Thus, P m in is 
not strictly the minimum of the dispersion relations but instead a point of symmetry - 
compare [126,91,69] and in particular section 5.4 below. 



The agreement with the prediction LP min = 2-7T is good for all values in table 3 bearing in 
mind that we have ignored systematic errors. This applies also to the point <j) = ip = ^p, 
A = | in the massless incommensurate phase because here the estimate for the error of L 
is probably too small. 



5.3. Correlation functions in the low-temperature regime 

In this section we apply the method explained in section 5.1 to the correlation functions 
in the low-temperature phase of the Z3-chiral Potts quantum chain. Note that the duality 
argument of section 3.2 applies only to the Hamiltonian and not to other operators. Thus, 
quantities like e.g. correlation lengths may be different in these two phases. 

We study again the correlation functions Cr(x) and C a (x) as defined by (5.1.2) and (5.1.5) 
using a low-temperature expansion for the ground state from the state | GS) (3.2.2). The 
expansion of the ground state in powers of A leads to an expansion of the correlation 
functions in powers of A: 

oo 

Cs(x) = ^A fc ci fc) (£). (5.3.1) 

k=0 

Below, we will first give the final results for general angles <p, p. We then specialize 
to the superintegrable case = ip = ^ and calculate even higher orders. By looking for 
a good fit we try to guess the structure of the correlation functions - as before. With 
this experience we turn back to the general case and discuss how the correlation functions 
should change for general cf>, p. 

In order to save space we present only the final results for the correlation functions. 
For Cr(x) one obtains, using the abbreviations C = cos (§),£= 1- AC 2 : 



Cf{x) =1 , 

cf\x) =-±-{^0-1} , 



6C 2 



4 1} (x) = 0, 
~ (3 ) cosy? 



1 



27C 2 



(1 - S x , ) 



2(1 - 16C 2 ) 1 



1 + 2C 2 - 3 i sin 



m 2 



16C 2 



+ 5. 



x.l 



(¥) 



m 2 



+ 



16C 2 



(5.3.2) 
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The first orders of the correlation function C a (x) read as follows: 

C<°>(aO=*.,o , CP(x) = 0, 

° K ' 9{C 2 ' \K 2C 2 J J CT 1 J 9C \2C 2 ft J 

^ 4) (-) = 

1 f / 7 8\ c ( 15 29 16C 4 \ r / 1 7 + 44C 2 



81C 2 [ V16C 2 Tl) ' \16C 2 12ft TZ 3 J ' \ 8C 2 12ft 2 

(5.3.3) 

Note that the correlation functions (5.3.2) and (5.3.3) do not depend on the charge sector. 

Eq. (5.3.2) suggests that the correlation function Cr(x) tends to a non-zero constant 
for large distances x - in contrast to the correlation functions in the high-temperature 
phase so that that the low-temperature phase is ordered over long ranges. However, beyond 
this general conclusion, it is difficult to guess from (5.3.2) or (5.3.3) what might be the 
behaviour even for small x. Thus, we set = (p = calculate four further orders and 
obtain 



4° 
4 2 



4 4 



4 6 



(x) =1 , (x) = c£ 3) (x) = (x) = 6jp (x) = , 

2 

fa) =g (*X,0 - 1) , 

0*0 = ^7 (^.o - !) + <W- 



81 v ' 7 ' 162 ' 
(x) =— {336 (4,o - 1) + 1604,1 + 604, 2 } - {264,i + 204, 2 } , 

4 8) 0*0 =^7^ i 12600 (<^° - 1) + 68525 X;1 + 35214,2 + 12254, 3 } 

- l ^§^ {9604,i + 9954,2 + 5254,3} 

and 

C<°> (s) =4,o , (*) = C< 3 > (*) = C< 5 > (x) = C<p (x) = 0, 

^ 2) W=^ {"44,0 + 4,i} , 

C?\x) ~ {-414,o + 144,i + 84,2} , (5.3.5) 

Cf\x) = {"5864,0 + 1474,i + 1264,2 + 804,s} , 

Cf\x) = ^777 {-99274,o + 21304,1 + 17214,2 + 17284,3 + 9104,4} • 

Up to the order calculated, C„(x) is real for all values of the parameters 0, <p, A. Cp(x), in 
contrast, has a non- vanishing imaginary part. By analogy to the high-temperature regime 
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and from the results in [105, 66] we expect that Cr(x) is oscillating. Indeed, the correlation 
functions in the superintegrable case = cp = ^ can nicely be fitted by 



C a {x) = ad xfi + be £ , 



Cr(x) = m 2 + pe 



■(£+ a e i )< 



(5.3.6a) 
(5.3.6b) 



For A G {j, |, | } good fits to (5.3.4) and (5.3.5) using (5.3.6) are given by the values in 
table 4. 



A 


6x 


a 


b 




m 2 


L 


V 


tnum. 




0.25 


0.26(1) 


0.891(9) 


0.099(9) 


0.24(3) 


0.9857605 


30 ± 18 


0.011(3) 


0.26(3) 


27 ±7 


0.50 


0.41(2) 


0.841(7) 


0.118(7) 


0.37(5) 


0.9381 


30 ± 15 


0.05(1) 


0.42(6) 


30 ±8 


0.75 


0.58(5) 


0.74(1) 


0.16(1) 


0.50(7) 


0.832 


30 ± 13 


0.15(2) 


0.6(1) 


28 ±4 



Table 4: Parameters for the correlation functions (5.3.6) at = ip = 

First, we remark that the correlation lengths satisfy £ CT = £r =: £ for all values of A within 
the numerical accuracy. In fact, one expects this equality because the correlation lengths 
should be the inverses of some mass scale, and there is only one mass scale in our problem 
because all three charge sectors are degenerate. Furthermore, we observe that our data 
is compatible with an oscillating correlation function for the operator V. The oscillation 
length L (or wave vector) is around 30 sites in a major part of the low-temperature phase. 
In [105, 66] it has been predicted that L should diverge as A crosses the phase boundary 
and approaches A = 1 where the critical exponent is expected to equal |. Our results are 
compatible with a divergent oscillation length at A = 1 although due to the large errors 
we do not even see that L increases with A. 

In addition to the perturbation expansions we can also perform numerical evaluations of 
the correlation function Cr(x). Since we can use vector iteration for up to N = 13 sites 
implying a distance of up to x = 6 one might hope to obtain more accurate results with 
this method. That this is not really the case is demonstrated by the values £f um ' and 
L num - which we have included into table 4. The main problem with numerical calculations 
is the numerical accuracy: In particular, the evaluation of the constant contribution m 2 is 
not free of errors. On the other hand the short correlation lengths strongly damp the non- 
constant part of the correlation functions: For example Cr(x) — m 2 has already decreased 
by 10 orders of magnitude at x = 6 for A = \. Therefore, for our estimates we have used 
only x < 4 for A = \ and x < 5 for A = \ and A = |. 

Thus, it will be very difficult to obtain more precise results from approximative arguments 
and an exact expression for Cr(x) is probably needed in order to decide whether (5.3.6b) 
really is the correct form and to determine the wave vector L accurately. 

Before we conclude the discussion of the correlation functions for the superintegrable chiral 
Potts model, we mention that a conjecture for the form of Cr(x) has been formulated in [4]: 
Cr(x) = m 2 + 0(e *F ) where m is the order parameter. Our result (5.3.6) is compatible 
with this from. In [97] the conjecture for the order parameter 



m = 



(v 


I 


X 


v) 


(v 


v) 



(5.3.7) 
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has been formulated, but (5.3.7) has not been proven yet. The constant term in (5.3.4) is 
in exact agreement with (5.3.6) and (5.3.7) up to the order calculated, such that we may 
assume that at least the constant term of Cr(x) is now known exactly. 

A few remarks on the choice of ground state (3.2.2) are in place because in [97, 11] 
(5.3.7) has actually been derived considering an expectation value of the operator r x . We 
have already pointed out that the one point functions of T x vanish identically due to charge 
conservation if one uses the charge eigenstates (3.2.2). However, if one uses instead non- 
charge eigenstates like | ... 0) for a perturbative expansion of | v) they do not vanish. 
Indeed, using an expansion for | v) from | ... 0) we once again verified equality of this 
one point function with the order parameter m. If we redefine Cr{x) by replacing \v) by 
| v) and subtracting the contribution from the one point functions, this is in fact the only 
change, i.e. Cr(x) = Cr(x) — m 2 . C a (x) remains unchanged under this redefinition. 

For more general values of the angles 0, ip one expects the correlation functions to be 
also of the form (5.3.6) - of course with different values of the parameters. We can see 
from the constant term m 2 of the correlation function Cr{x) (5.3.2) that it will not be of 
the form (5.3.7) for general ^ ^ ^ cp. In general, the coefficient of A 3 for the constant 
term does not vanish and m 2 does not even have an expansion in powers of A 2 . Among 
the powers that we have calculated for the general case only the fourth order in (5.3.2) 
has a non-vanishing imaginary part at x = 1. Under the assumption that (5.3.6b) is the 
general form we would expect the imaginary part at x = 1 to be proportional to sin (^) 
for very small temperatures A. Thus, we expect for very low temperatures A the relation 
L _1 ~ (p. On the one hand, this explains the conjectured presence of a second length scale 
L in addition to the correlation length The oscillation length L just comes from the 
chiral angle 4> and thus these two scales must be related to each other. On the other hand, 
the oscillation (should it really be present) will vanish smoothly as the parity conserving 
Potts case = (p = is approached. 



5.4. Form factors and correlation functions 



In this section we show how to use form factors in order to determine the oscillation 
length L exactly using symmetries of the Hamiltonian arising for particular values of the 
parameters. 

One can use the quasiparticle picture which we have discussed earlier in order to rewrite 
a correlation function Cz(x) as follows: 



C s (x) 



E~= jiT (nr=i d Pi) ( v 


■=■+ 


|pi,...,Pn) (pi,...,Pn| 


Si 


v) \(v\ 




h 


v)\ 2 


(v\ 


V) 






M 


v) 


2 





n=l" u \i=l 



(5.4.1) 

where we have inserted a complete set of normalized n-particle states \p±, . . . ,p n ) . Rep- 
resentations similar to (5.4.1) have been used in quantum field theory for a long time (see 
e.g. [107]) and are well-known to be useful for the evaluation of correlation functions of 
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statistical models (see e.g. [111]). According to (5.4.1) one could compute the correlation 
function C~,(x) by computing its L form factors' (p±, . . . ,p n | Si | v) , but one can even 
derive interesting results without doing so. Clearly, if the groundstate | v) has non-zero 
momentum ^ we expect an oscillatory contribution to the correlation function. 
However, one can read off from (5.4.1) that an oscillatory contribution is also to be ex- 
pected if P^ = but the model breaks parity which precisely applies to the massive 
high-temperature phase of the chiral Potts model. The correlation functions of massive 
models in general have an exponential decay, i.e. Cs(x) = e~ ^ f~.(x) where f~.(x) is some 
bounded function. According to (5.4.1) we also expect an oscillatory contribution of the 
form e % ~T^ . In summary, we expect correlation functions of the approximate form 



271-3: 



C s (x) ~ e -« +l — . (5.4.2) 

£ is called 'correlation length' and L is the 'oscillation length' (L -1 is the 'wave vector'). 
Eq. (5.4.2) explains the fits (5.1.10) and (5.3.6) that we have used earlier. More precisely, 
for the Z n -chiral Potts model the operator Ti creates Q = 1-single-particle excitations 
from the groundstate. The dispersion relations of these particles clearly violate parity. 
Therefore we expect that Cr(x) is of the form (5.4.2). The action of the operator o~i is 
much less spectacular. In particular, it leaves the charge sector Q = invariant and thus 
it need not necessarily have an oscillatory contribution. In fact, from (5.1.3) we see that 
C a (x) should be real which in view of (5.4.2) implies the absence of oscillations. 

For massive perturbations of conformal field theories it is well-known that form factor 
expansions already for few particles yield excellent approximations (see e.g. [111]). In 
order to gain some intuition we perform some simple computations for the form factors of 
the operator V for the Z3-chain. 

First, considering perturbation series for the states one sees that 

(p u ...,p r \ r x \v) =0(\M). (5.4.3) 

Next we compute the states \v) and \ p±) from perturbation series of (4.2.1) and (4.2.2) 
respectively. One finds that up to order A the states are already normalized properly, i.e. 

(v\v) = 1 + 0(A 2 ), (pilpi) =1 + 0(X 2 ). (5.4.4) 

Finally one finds the single-particle form factor for n = 3 up to order A: 



(5.4.5) 



Comparing (5.4.5) with (4.2.7) up to order A one sees that the maximum of the form 
factor is obtained when Ei(pi) has a minimum and vice versa. Furthermore, multi-particle 
states have higher energy than single-particle states and due to (5.4.3) the contributions 
of the form factors decrease with the number of particles for small A. In summary, one 
obtains the qualitative statement that the contributions of the form factors decrease with 
increasing energy. 



i cos \pi — ■& 

(Pil Fiji;) =- j= ll + X V j— I +0(A 2 ) 
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According to (5.4.3) the single-particle form factor (5.4.5) is the only one contributing to 
C-p(x) up to order A. Indeed upon inserting (5.4.5) into (5.4.1) one recovers (5.1.7) up to 
first order in A - as it should be. Note however that the direct computation of (5.1.7) is 
simpler than the one using a form factor expansion. 

Symmetries of the Hamiltonian translate into symmetries of the form factors. In 
certain cases these symmetries are already sufficient to compute the oscillation length L 
as we will see now. 

Def. VI: (Parity operator) The parity operator *p is defined by the following action on the 
states (3.1.6): 

r(ty) \h in) =\hiNiN-i h) ■ (5.4.6) 

Remark: ^cri +x *P = cri_ x , ^STi^ty = Ti- X . 

The following lemma is based on results obtained via numerical diagonalization in [126] 
which showed that the Hamiltonian (3.1.1) has exact symmetries for particular values of 
the parameters: 

Lemma I: Denote the restriction of the Hamiltonian Hffl in eq. (3.1.1) to the spaces 

with momentum P and charge Q by 'H^\p, Q) '. Then one has the following identities 
(see also [126]): 

a k = a n - k ^H^\p,Q)^ = H^\-P,Q), 

a* k = a n - k anda k eR ^H^ ] (P, Q)<$ = (h™ (-P, Q)) + , 

a* k =a n - k andal = e- 2 ^ zk a k yH$\P m , Q + P, Q)<$ = (H$\P m , Q - P, Q)) + 

(5.4.7) 

where the symmetry of the last line holds for any Q that is invertible in Z n and those P m ,Q 
satisfying -P m ,Q<5 _1 + ttz = mod tt as well as e i2Pm - Q being an nth root of unity. For very 
particular values of the momentum P and the chain length N some projection is necessary 
in (5.4.7). m 

The first two lines of (5.4.7) follow immediately by looking at tyH^ty, keeping in mind 
that the translation operator defined in (3.1.10) satisfies tyT N ty = T^ 1 = T^. The 
derivation of the third line of (5.4.7) is more complicated and relies on a suitable choice of 
basis. Because of its technicality we have shifted it to appendix B. 

Remarks: 

1) For the last line of (5.4.7) insert the parameterization (3.1.12) and set z = Then 
Pm,Q = tt(1 — ^7) is a solution to -P m ,Q<5 _1 + = lying in the interval [— tt, it] 
the other solution is shifted by tt. The solution P m ,Q = 7r(l — — ) corresponds to the 
minimum (4.2.4) of the dispersion relation of the single-particle state in this charge 
sector. 

2) The restriction to Q invertible in Z n (those Q whose greatest common divisor with 
n is 1) was necessary for technical reasons but there are indications that lemma I is 
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true without this restriction. Combining this with remark 1) we conclude that the 
relation P m ,Q = tt(1 — — ) seems to be valid for general Q. In fact, close inspection 
of the argument in appendix B shows that the assumptions can be weakened in order 
to cover e.g. also the case Q = 2 for n = 4. 
3) In the generic case (which is the only case we are really interested in) one may ignore 
the projection operators in the equality (5.4.7). 

Theorem VI: (Oscillation length [91]) Let the Hamiltonian H(P, Q) restricted to momen- 
tum and charge eigenspaces with eigenvalues P and Q have one of the following symmetries: 

yH(P m , Q +P, Q)<p = H(P m , Q -P, Q) or ^H(P m , Q +P, Q)% = (H(P m , Q - P, Q)) + 

(5.4.8) 

with some P m ,Q depending on the charge sector Q. Assume furthermore that P^ = and 
that Si \v) has charge Q. Then the oscillation length L of the correlation function C~,(x) 
satisfies 

LP m , Q = 2tv. (5.4.92 

Proof: We start from the form factor expansion (5.4.1) which in the present case becomes 

C*{x) = Y &Pe iPx l(P,Q; / r| , Ellv)l " (5.4.10) 
^Jo (v\v) 

where we have only written the quantum numbers P and Q explicitly and comprised the 
other ones in the label 'r'. First we observe that ^3£i*J} = Hi. If the Hamiltonian satisfies 
< ^H(P ult Q+P, Q)^P = H(P m> Q—P, Q), then eigenstates of momentum P m ,Q+P are mapped 
under parity to eigenstates of momentum P m ,Q—P- This means that ((Pm,Q+P), Q',t\ Hi | 
v) = ((P m ,Q — P)i Q'i f | Si | v) . If the symmetry involves the adjoint of the Hamiltonian 
one finds ((P m ,Q + P),Q;r\ Hi | v) = ((P m ,Q ~ P),Q;r\ Hi \ v)*. Thus, the following 
identity is valid in both cases: 

| ((P m , Q + P), Q; r | Hi | v) | 2 = \((P m , Q - P), Q;r\ Hi \v)\ 2 . (5.4.11) 

Now we turn back to the form factor expansion (5.4.10): 

C s (x) = y < J dPe^ mQ; ^ llv)l2 + I dPe 1 

Pm,Q Pm,Q— 7T 



\(P,Q;r\ 


Si 


v)\ 2 


(v 


\v) 



\{(P m , Q + P),Q;r\ 


"1 


v)\ 2 


(v\ 


v) 



+e 



\((P m , Q -P),Q;r\ 


= 1 


v)\ 
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(v\ 


v) 
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e iPm, Q xJ2 J dP2cos(Px) 
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(v\ 


v) 



(5.4.12) 
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where the last equality follows from (5.4.11). This shows that C~,(x) is of the form 



C s {x) = e^f{x) (5.4.13) 

with L satisfying (5.4.9) and f(x) is given by the remaining integral in (5.4.12) which is 
clearly real. □ 

Remarks: 

1) Note that theorem VI is true for more general Hamiltonians H(P, Q), but it covers in 
particular the case (5.4.7) for the Z n -chiral Potts model. 

2) Note that if the Hamiltonian has several different -P m) Q such that (5.4.8) holds (which 
applies to (5.4.7)) one obtains different expressions for C'e(x) involving different L 
and f(x). We have argued before that form factors are maximal for smallest energy. 
Therefore one may expect to obtain the representation with f(x) > on x G [0, y] 
with maximal y choosing the value P m ,Q that corresponds to the minimum of the 
dispersion relation of the fundamental particle. 

From lemma I and theorem VI we immediately obtain the following 

Corollary: 



C rQ (x) =e — f Q , r (x) (5.4.14a) 



with 



f Q , r (x)eR Vx, 

L = oo for (p = tti, r E Z or (p e R, Re(^) = 



In 



(5.4.14b) 



n-2Q 



for pel, Re(0) = 7rand0<(5<n 
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// libro del la natura e scritto in lingua matematica 

Galileo Galilei 



6. Extended confer mal algebras 

Symmetries play a fundamental role in theoretical physics because they help to simplify a 
problem considerably or sometimes they even determine the relevant quantities completely. 
For example, the symmetry group SO (4) of the quantum mechanical central force problem 
allows to determine the spectrum of the hydrogen atom by purely algebraic computations. 
In particular, for quantum field theories it is desirable to have a sufficiently large symmetry 
because they can be solved directly in only very few cases. Such a situation does indeed 
arise in two space-time dimensions where so-called 'extended conformal algebras' (also 
called W-algebras) provide infinite symmetries. The study of the representation theory 
of these algebras gives rise to 'rational conformal field theories' (RCFTs) which are in 
some sense finite. Particular RCFTs describe the second order phase transitions in Z n - 
spin chains as we have already mentioned earlier. More generally, a classification of RCFTs 
would also classify the associated conformally invariant critical points of statistical models. 
Therefore, it is interesting to understand the structure of RCFTs. This also involves a good 
understanding of W-algebras and the RCFTs associated to them (their so-called 'rational 
models'). The remainder of this thesis is devoted to investigations that either yield new 
results on the structure of W-algebras or are relevant for spin quantum chains. 

6.1. Quantum W-algebras 

In this section we give a brief review of the notion of a quantum W-algebra - for details 
see e.g. [33, 55]. 

The ingredients of a quantum W-algebra are: 

1) A quantum W algebra is a vector space of local chiral quantum fields, i.e. there is a 
basis of fields <f>(z) depending on a complex coordinate z G C (including one point at 
infinity) with definite 'conformal dimension' d(4>) and mode expansion 

Cj>{z)= ^'^^n (6.1.1) 

n-d(»eZ 

such that all fields <p(z) are either bosonic or fermionic and that their dimensions are 
half-integer 2d(4>) G Z. 

2) The algebraic structure is encoded in the so-called 'operator product expansion' (OPE) 
of two fields (f>(z), x( w ) : 

ib(w) 

^)°xM = E ( ,_ ro) S +J M-, W (6-1-2) 

where the sum on the r.h.s. of (6.1.2) runs over all local chiral fields in the W-algebra. 
The A^ x are the structure constants of the W-algebra. The l.h.s. of (6.1.2) is assumed 
to be radially ordered. The OPE (6.1.2) is required to be associative. Inserting the 
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mode expansion (6.1.1) into (6.1.2) and performing Cauchy integrals, the singular 
part of the OPE gives rise to a Lie bracket whereas the regular part defines a normal 
ordered product ((pd n x) (the derivative d is defined in the obvious manner). Note 
that the associativity of the OPE can be used in order to compute OPEs involving 
normal ordered products (also known as 'Wick rules'). 

3) The W-algebra contains a distinguished field, the energy-momentum tensor L with 
OPE: 

. * . . c 2L(w) L'(w) 

L(z) o L(w) = ^ + 1 2 + r -^- + reg. 6.1.3 

2{z — (z — w) z (z — w) 

The constant c is called the 'central charge'. The Lie algebra obtained by a mode 
expansion of the singular part of (6.1.3) is a representation of the 'Virasoro algebra' 
where the central element C is represented by cl. 

4) The energy-momentum tensor L defines a grading on the W-algebra: 

d(<j>)<p(w) ^ ^ Aj^jw) 



where d(cp) is the conformal dimension of 
5) Fields transforming as 



L z o to = — — \- - + y - ~7/ 1 \ (6.1.5) 

W V 7 (Z-W) 2 (Z-W) ^ (2 - u;)2+d(0)-d(V) v 



|d(V)-d(*)|>l 



are called 'quasi-primary' . The field is called 'primary' if the remaining sum runs 
only over fields ip with d(ift) > d((p) + 1. Note that the energy-momentum tensor L 
itself is quasi-primary with conformal dimension d(L) = 2. 

6) A basis for all fields in the W-algebra can be obtained from a countable basic set of 
quasi-primary fields {W^} using normal ordered products and derivatives. The 

are called 'generators' . If the set of generators is minimal, the W-algebra is called a 
l W(d(wW),d(wW),...)\ 

7) There is an irreducible representation of the mode algebra, the so-called 'vacuum 
representation' , i.e. the modes 4> n of any field <f>(z) act as endomorphisms in the vacuum 
representation. The vacuum representation has a cyclic highest weight vector | v) with 
the property 

cf) n \ v ) = 0, Vn < d(<f>) . (6.1.6) 

A normal ordered differential polynomial that acts trivially in the vacuum representa- 
tion is identically equal to zero. Fields with this property are called 'null fields'. Note 
that it is sufficient to require (6.1.6) for quasi-primary fields only, it then automatically 
extends to all fields including derivatives. 

Property 7) gives a precise meaning to the fields <p{z). The definition of a W-algebra can 
be made precise using the mathematical language of vertex operator algebras (see [47] for 
a definition). Then, property 7) is the starting point of the definition and properties 1) - 
6) turn out to be reformulations of the axioms of a vertex operator algebra. In general, one 
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will have to impose additional finiteness conditions like requiring the graded components 
of the vacuum module to be finite dimensional. 



The simplest W-algebras are 'current algebras' '. They are generated by currents Ji(z) 
(fields with dimension d{ Ji) = 1) with OPEs 

J i (z)oJ j (w) = - — ^ + ) f V +reg. 6.1.7) 

(z — w) z ^ (z — W) 

where gij is an invariant metric of some Lie algebra C n and the are the structure 

constants of this Lie algebra. The current algebra defined by (6.1.7) is denoted by '£ n ' 
and k is called the L level'. The singular part of the OPE of the basic currents Ji defines 
an infinite dimensional Lie algebra which is a so-called 'Kac- Moody algebra'. Below we 
will abuse notation and also refer to the W-algebras defined by (6.1.7) with the term Kac- 
Moody algebra. In particular we will ignore that it is not entirely trivial that a Kac-Moody 
algebra becomes a W-algebra once it is endowed with a normal ordered product. 

There are several different constructions for W-algebras, most of them starting from 
Kac-Moody algebras (6.1.7). The so-called l coset construction' (also known as the GKO 

construction since for t n = sl(2) it is the construction for representations of the Virasoro 
algebra (6.1.3) of ref. [78]) is particularly simple to formulate. There are two closely related 
types of cosets: 

W W 

(6-1.8) 



n 



By definition, the coset W-algebra Wj C n is formed by all fields in the W-algebra W that 
commute (i.e. have regular OPEs) with the Kac-Moody subalgebra C n C W. Similarly, 
the coset W/C n is defined as the set of all fields in W that are invariant with respect to 
the action of the Lie algebra C n implemented by the zero modes of the currents t n C W. 
If the energy-momentum tensor does not coincide with the Sugawara energy-momentum 
tensor of £ n , the coset (6.1.8) yields automatically a closed W-subalgebra of W (see e.g. 
[33] for more details). 

Another very important construction is the so-called 'quantized Drinfeld-Sokolov' (DS) 
reduction (see e.g. [32, 60] for recent descriptions). This reduction is a Hamiltonian reduc- 
tion that also gives a precise meaning to the 'quantized Miura map' of the Fateev-Luk'yanov 
construction (see e.g. [54, 55]). Here we just recall some ingredients - for details see e.g. 
[32,60]. The data specifying the DS reduction is an embedding of si (2) into C n . The 
Cartan generator H of the sl(2) defines a grading on C n . By definition, the positive step 
operator E + of this sl(2) has grade one. Let T J be the basis of C n with the metric gij 
and structure constants /■■ used in (6.1.7). We introduce a Lie algebra valued current 
(J) by (J) = J2j JjT 3 ' ■ The grading on C n leads to the following decomposition of the 
Kac-Moody algebra 

(j) = (j)_ e (j) e (J)+ (6.1.9) 

where ( J)o refers to the currents of grade zero, (J)- to all currents with arbitrary negative 
grade and ( J)+ to the currents with arbitrary positive grade. The DS reduction is defined 
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by imposing the following first class constraints on the currents with positive grade 5 ): 



(J)+ = E+. (6.1.10) 

This amounts to equating the positive grade currents to zero except for the one corre- 
sponding to E + that is set to one. From a classical point of view the constraints (6.1.10) 
generate gauge transformations and the problem is to find a suitable representative for each 
gauge orbit. At the quantum level the constraints (6.1.10) are imposed using a nilpotent 
BRST operator Qbrst (Qbrst = 0)- I* 1 order to impose the constraints with a BRST 
operator one has to introduce some auxiliary fields, so-called ghost-antighost pairs. Then 
one considers the cohomology of Qbrst and the problem is to find a suitable representative 
for each cohomology class. 

Suitable representatives are associated to the so-called lowest weight gauge: The generators 
of the DS W-algebra can be chosen to be the lowest weight components of the 
embedded sl(2) in the adjoint representation on the Kac-Moody algebra C n . From the 
lowest weight gauge one reads off that the conformal dimension of W^' is given by the 
spin Sj of the si (2) representation, i.e. d(W^) = Sj + 1. Furthermore, the generators can 
be made primary using a certain 'dressing'. In the quantum case, this dressing is provided 
by the 'tic-tac-toe' equations of [32]. 

For the principal embedding of sl(2) into C n the DS reduction gives rise to the so-called 
Casimir W-algebras which we denote by W£ n . In the case of a nonprincipal si (2) embed- 
ding into a simple Lie algebra C n we use the notation W§ n where S denotes the embedding. 
For example, for C n = A n , S can be chosen as the r-tuple of the dimensions of the ir- 
reducible sl{2) representations which appear in the defining representation of C n . The 
Polyakov-Bershadsky algebra W(l, §, |, 2) [116, 25] which is obtained by DS reduction of 



the nonprincipal sl(2) embedding into si (3) is abbreviated in this notation by W 



sl(3) 
2,1 • 



In the next section we shall be interested in non-trivial outer automorphisms of W- 
algebras and their effect on the HWRs. An automorphism p of a W-algebra is a bijective 
map of the algebra that satisfies the following compatibility condition with (6.1.2): 

,(«*)) o pfa W) = Y, (z -^wtS-dw • < 6 - L11 > 

p is called an 'outer' automorphism if it is not generated by the W-algebra itself. Each 
such automorphism enables one to impose non-trivial boundary conditions on the fields <f)j 
in the algebra: 

4> j (e™z)=p{<l> i {z)) (6.1.12) 

where now the argument z lives on a suitable covering of the complex plane. This type of 
boundary condition will be called a 'twist'. 



5 ) There are some subtleties for the grade \ part of a half-integral embedding which we 
ignore here because it is not relevant for the aspects we need below. 



52 



Quite often automorphisms p have the form p(Wj(z)) = ±Wj(z) where some genera- 
tors Wj(z) transform with a plus sign and others with a minus sign. For any field with 
4>j (e 27 "z) = —(f) j(z) the Laurent expansion (6.1.1) has to be modified to 

M*)= E * n_d(< ^>. (6-1.13) 

n-d(^)eZ+i 

This leads to the Ramond-sector of a fermionic W-algebra. Bosonic fields acquire half- 
integer modes. 

In the presence of an outer automorphism one can also consider the projection of the 
W-algebra onto the invariant subspace. This projection is called an ' orbifold . 

6.2. Automorphisms of Casimir W-algebras and boundary conditions 

The simplest W-algebras are those with one additional generator, so-called ' W(2, 8)- 
algebras'. The automorphisms of these algebras were studied in detail in [90] and [88]. 
Here we will focus on the so-called 'Casimir' algebras W£ n where the dimensions of the 
simple fields equal the orders of the Casimir invariants of a simple Lie algebra C n (see 
e.g. [15,16]). One common approach to their study is Toda field theory [19,20,26]. As 
we will see below only Casimir algebras W£ n based on simply-laced C n have non-trivial 
outer automorphisms for generic c. The unitary minimal series of these algebras can also 
be studied via GKO-constructions [78]. Q. Ho-Kim and H.B. Zheng have noticed that in 
this approach outer automorphisms of the Lie algebra give rise to automorphisms of the 
W-algebra and argued that there no further ones [85,86,87]. Owing to their work twists 
of the unitary minimal series of Casimir algebras are well understood. Still, we would like 
to comment on Casimir algebras from the point of view of extended conformal algebras, 
including the non-unitary minimal models into our discussion. 

One way to study automorphisms of Casimir W-algebras is by inspection of the struc- 
ture constants of known examples. Inspecting the structure constants of the algebras 
constructed e.g. in [30, 102, 93, 94] it was observed in [90] that W£ n with generic central 
charge c has as many automorphisms as C n has. This was also observed before for the 
unitary minimal series in [85,86,87]. In fact, it can easily be seen that using quantized 
Drinfeld Sokolov (DS) reduction (see e.g. [32,60]) each outer automorphism of the un- 
derlying Lie algebra gives rise to an automorphism of the associated Casimir W-algebra. 
In the DS framework Casimir W-algebras arise from the principal sl(2) embedding. The 
Cartan- resp. step-operators of this principal s/(2) are sums over all the Cartan- resp. 
step-operators associated to the simple roots of the underlying Lie algebra. Since outer 
automorphisms of simple Lie algebras are nothing but certain permutations of the roots 
the principal sl(2) is clearly preserved. Therefore, the automorphism of the underlying Lie 
(or Kac-Moody) algebra survives the reduction. We will argue below for WA n that the 
automorphism remains non-trivial after the reduction. In the case of Casimir W-algebras 
it seems that this type of automorphism is in fact the only one [90] although this is far 
from being obvious. In fact, also W-algebras for non-principal s/(2) embedding admit 
automorphisms some of which do not seem to have any relation to the automorphisms 
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of the underlying Lie algebra (compare e.g. the analogue of the Ramond sector of the 
Polyakov-Bershadsky-algebra W(l, f , f , 2) [116, 25]). 

In the case of WA n the effect of the automorphism of the underlying Lie algebra after 
reduction can easily be controlled. To this end one uses the lowest weight gauge for the DS 
reduction, i.e. one has to find the lowest weights of the principal sl(2) in the adjoint repre- 
sentation. These lowest weights are given by the tensor product of a spin ^ representation 
with itself minus a spin zero representations (in order to ensure tracelessness). By inspec- 
tion of the Clebsch-Gordan coefficients one sees that the parity under the automorphism 
alternates with the spin S, i.e. the lowest weights transform with a factor (— l) s . Since 
the dimension of the generating fields W^' is d(W^) = Sj + 1 and the energy-momentum 
tensor L = must be invariant, one sees that the generators with even dimension are 
invariant and those with odd dimension pick up a minus sign. This shows that 

P (WW) = (-l) d(w, "V (j) (6.2.1) 
for the generators of WA n = W(2, . . . , n + 1). 

For certain values of the central charge one finds null fields in the algebras WA 2 — W(2, 3) 
and WA3 = W(2, 3, 4). These can be exploited to determine the conformal dimensions of 
some minimal models along the lines of [48, 123] using the special purpose computer algebra 
system [89]. Certain normal ordered products in the twisted sector are not straightforward 
to handle (see e.g. [50] for a discussion of problems of this type) only a few models are 
accessible to these direct computations. The results have been presented in [90]. Here, we 
list only the subset of results corresponding to Z n parafermions (see table 5). 



W(2,3) 


<= = ! 


W(2,3,4) 


C= 1 


untwisted 


twisted 


untwisted 


twisted 





1 





3 


1 


8 


4 


16 


2 


1 


1 


9 


3 


3 


40 




16 


16 


2 




1 


1 


1 


5 




3 


16 


48 


1 




1 




25 


15 




12 




48 



Table 5: Conformal dimensions of untwisted and twisted 

W(2, 3) = WA 2 and W(2, 3, 4) 9* WA 3 for the Z 3 resp. 
Z 4 parafermionic models c^ 2 (4,5) = |, c^ 3 (5,6) = 1. 

From the above discussion we conclude that there will be no outer automorphisms for 
WAi, WB n: WC n , W£ 7 , W£ 8 , and WG 2 . The algebras WA n for n > 1, WV n 

for n > 4 and W£q should have exactly one outer automorphism and correspondingly 
exactly one twisted sector in addition to the untwisted one. For these algebras explicit 
computations [90] showed how the formulae for the h- values in [86, 87] might generalize: 

cc n (p, q) 

/-(l, 2, 3) 

h " 



= n- Yip- 



(p - qY 



pq 



(pX - qp) : 
2pq 



+ 



r 



n 



(6.2.2) 
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where 



n no no 

i=l i=l i=l 



(6.2.3) 



and no, n\ are the dimensions of the invariant subalgebra C' no respectively twisted subal- 

gebra of C n ; Qi are the fundamental weights of C n ; Qi the fundamental weights of C' n() ; h 
is a constant (that can be interpreted as the conformal dimension of some field) and r^, Si 
positive integers subject to certain constraints. For the case of A n the invariant subalgebra 
is Cp+i j. In the unitary minimal series of WA n , WD n (n > 4) and WSq one has h = j^. 

We observe that (6.2.2) indeed reproduces the data in table 5 for W(2, 3, 4) if we use the 
weights of C 2 (which is the invariant subalgebra) and h = jq. 

For C n (periodic boundary conditions) eqs. (6.2.2) with (6.2.3) have been rigorously 
proven in [63] (see e.g. [28] for a brief summary including the range of the integers r« and 
Si). It should be possible to generalize the result of [63] to twisted Kac-Moody algebras 
Cn^ and yielding a rigorous derivation of (6.2.2) with (6.2.3). At least one thing is 
immediately clear if one assumes that the representations of a (twisted) W-algebra obtained 
by DS reduction are induced by admissible representations of the (twisted) Kac-Moody 
algebra: The representations of a Kac-Moody algebra are labeled by the weights of the 
Lie algebra formed by the zero modes of all currents. While for an untwisted Kac-Moody 
algebra Cn^ the algebra of zero modes is isomorphic to C n , the zero modes of a twisted 
Kac-Moody algebra form a Lie algebra C' nQ that is precisely the invariant subalgebra of 
r 

The exceptional cases W£§ and WV4 are particularly interesting. Especially for 
WT>4 = W(2,4,4, 6) the group of outer automorphisms should be S3 for generic c. In 
W(2, 4, 4, 6) one structure constant remains a free parameter [101]. It was observed in [90] 
that one can define an operation of 0(2) on the two generators of conformal dimension 4 
under which the structure constants transform covariant. The algebra is invariant under 
the natural embedding of S3 into the group 0(2) operating on the fields if one chooses the 
self coupling constants of these two generators to be equal. 

Denote the primary fields of dimension 4 in W(2,4,4, 6) by V(z) and W(z). Then the 
iS3-symmetry of this algebra translates into the following type of boundary conditions for 
the given choice of coupling constants: 



to those elements of S3 which under the embedding yield elements of 50(2). The boundary 
conditions (6.2.5) correspond to the three elements of S3 that are mapped to elements in 
(9(2) with determinant —1. 



V(e 2 ^z) 
W{e 2 ^z) 



cos(a)V(z) — sin(a)W(z) 
sm(a)V(z) + cos(a)W(z) 



(6.2.4) 



or 
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The /i-values in the unitary minimal series (6.2.2) of WV4 have been calculated in [86] 
without having to consider the boundary conditions of the additional simple fields which 
look quite strange at first sight. Note that for V4 the invariant subalgebra is the exceptional 
algebra Q2 and here h in (6.2.2) satisfies h = jg. 

In [33] it has already been stated that the iS3-symmetry should lead to modes in -3 and the 
precise statement was presented in [90] which we now would like to recall. Let us first focus 
on (6.2.4). Set U^(z) := V(z) + iW(z) and U^(z) := V(z)-iW(z). Then (6.2.4) turns 
into U^\e 27Ti z) = e ia U^\z) and U {2 \e 2 ™z) = e" iQ t/ (2) (z) which can be satisfied by 
choosing modes in Z + 7p- for and those for in Z— Consider now (6.2.5). For 
this case set (z) := cos(a)V(z) + (sm(a)+l)W(z) and Y^(z) := cos(a)V(z) + (sm(a)- 
l)W(z). Now (6.2.5) turns into Y^\e 2ni z) = Y^\z) and Y^ 2 \e 2 ^z) = -Y^ 2 \z). This 
can be satisfied by choosing the modes for Y^ 1 ' in Z and those for Y^ in Z + |. 

We would like to conclude this discussion by commenting on the relation of twisted 
representations to orbifolds (for a detailed discussion of orbifolding see e.g. [44]). For W- 
algebras with a Z2 automorphism p 2 = 1 one has two partition functions, one where only 
the characters x W = tryq L °~^ of the untwisted sector enter, and one where the characters 
X W and x w = try pq L °~^ of both sectors enter. The latter can be identified with the 
partition function Z of the orbifold 6 ): 

z= E VxVrxV+(xVn?]+ E luvrxV+ffirxV] 

k: untwisted k: twisted 

o \ ^ ri/ w 1 ~w\*it w 1 ~W\ 1 1/ w ~w\*it w ~W\~\ (6.2.6) 

= 2 2^ 1.5 ufc +Xk) hbck +Xk) + i(xk -Xk) i(xk -xDJ • 

k: untwisted 
and twisted 

This implies that the characters of the orbifold W-algebra are given by §(XaT Xk") an< ^ 
^(xY — Xk)- We conclude that the h- values for the HWRs of the orbifold are those of 
the original W-algebra in both sectors in addition to some which differ by (half-) integers. 

It would be interesting to generalize the observations of this section to the supersym- 
metric case. Of special interest is osp{A\A) which is the supersymmetric analogon of V 4 . 
The corresponding Casimir algebra is a 5W(|, 2, 2, |). In [90] it was argued on the basis 
of some explicit results in [29] that 5W(|,2,2, |) should be invariant under the natural 
embedding of £3 into the group 0(2) operating on the fields if one chooses the self cou- 
pling constants of the two generators of dimension two to be equal - precisely like for 
yy£>4 = W(2,4,4,6). 

6.3. Consequences for spin quantum chains 

As already pointed out in [85] there is a close connection of the partition function of 
W(2, 3) at c = I including the twisted sector and the three states Potts model. In fact, the 
different boundary conditions of W(2, 3) correspond to the different boundary conditions 
of the Potts quantum spin chain at critical temperature. Choosing the spin shift operator 
at the end of the chain to be equal to the spin shift operator at the first site yields the 

6 ) We simplify notation by absorbing multiplicities of characters into the index set. 
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field content of the 'untwisted' sector of W(2, 3) (see e.g. [57, 133] and [71] for numerical 
verification). We recall the remarkable fact [85] that the twisted sector of W(2, 3) yields 
additional representations which can be identified with fields in the thermodynamic limit 
of the three states Potts quantum spin chain if the spin shift operator at the end of the 
chain is chosen to equal the adjoint of the one at the first site [37] . 

It would be interesting to know if this observation generalizes to all Z n . For twisted 
boundary conditions only partial results are available in the literature (see e.g. [74]). We 
shall therefore present an explicit verification of this statement in the case of Z 4 . We will 
follow the approach of [71] and study the spectrum Hamiltonian (3.1.1) at = ip = 
numerically. 

(n) 

Let i£jv,i be the eigenvalues of H N with periodic boundary conditions in ascending order 
and En,i those with twisted boundary conditions. Then the relevant scaling functions are 
given by [36, 72]: 

N ~ 

CJV ?' •= (Em j — Em n ) 

2tt ' ,u; (6.3.1) 

& : = J im Cam- 

In the case of periodic and cyclic boundary conditions, the eigenvalues Q of the charge 
operator (3.1.9) and momentum P are good quantum numbers. In the case of twisted 
boundary conditions neither charge nor momentum are conserved any more and one does 
not have any obvious conserved quantities. At least for even n the charge Q mod 2 is 
conserved. 



Assume that (3.1.1) with = ^ = exhibits conformal invariance at A = 1 and denote the 
dimensions of the fields in the left chiral part by h and of those in the right chiral part by 
h. For periodic and twisted boundary conditions the field theory is diagonal, i.e. the fields 
<f>(z, z) with dimension h + h satisfy h = h and thus have vanishing spin h — h. Therefore, 
the modes of the fields <f>(z, z) yield levels in the spectrum with £ = h + h + r where r G Z + 

for periodic boundary conditions and r G for twisted boundary conditions. 

In order to test this method we shall first study the well-known three states Potts model. 
For Z3 we have studied 3 to 9 sites. Thus, it was necessary to partially diagonalize matrices 
of dimension 3 9 = 19683. The limits — > 00 of the lowest gaps £j are given in table 6. 



z 3 


% 


ii 


h + h + r 





0.050000(2) 


— + — 

40 T 40 


1 


0.250005(5) 


8^8 


2 


0.5500(5) 


— + — + 1 

40 ^ 40 ^ 2 


3 


1.05(4) 


-J- + - + 1 
40 ^ 40 ^ ± 



Table 6: The low- lying spectrum of the 
twisted Z 3 -chain at A = 1 
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The numbers in brackets indicate the estimated error in the last given digit. For details on 
the extrapolation procedures and error estimation see e.g. [84] . We do not give more than 
four levels because the errors of the next levels make an accurate identification impossible. 
Note that we can nicely identify the dimensions ^ and g of the chiral fields - as expected. 

Let us now turn to Z4. The spectrum of the Z4-chain was already discussed in [14] 
also for twisted boundary condition. This was done applying Kac-Moody algebra and 
numerical techniques. Since the former considerations are analogous to the discussion of 
the previous section and the numerical results are unpublished, we will present results of a 
direct calculation here because we would like to demonstrate the correspondence between 
boundary conditions in statistical mechanics and conformal field theory. Note that the 
Z4-version of (3.1.1) at = <p = is a special case of the Ashkin- Teller quantum chain 
which was introduced in [98] setting the parameter h = | (in the notations of [13]). 

For Z 4 we have at least a splitting of the spectrum into two sectors of Q mod 2. We have 
studied 4 to 8 sites, implying the partial diagonalization of matrices in dimensions up to 
4 8 /2 = 32768. 
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h + h + r 


i 




h + h + r 



1 

2 


0.04167(2) 
0.375(2) 
1.040(2) 


48 T 48 
3,3 
16 16 
25 1 25 
48 48 



1 

2 
3 


0.1254(1) 
0.6231(3) 
0.623(7) 
1.12(1) 


16 T 16 

J_ + J_ + I 

16 T 16 ~ 2 

— + — + 1 

16 T 16 ~ 2 

^ + ^ + 1 

16 T 16 T 



Table 7: The low-lying spectrum of the twisted Z 4 -chain at A = 1 

The dimensions jg, jq and || of the chiral field theory can be nicely seen in the explicit 
results of table 7. 

The field content of the Z n -spin quantum chains (3.1.1) at their second order phase 
transition (j) = ip = 0, A = 1 is given by the first unitary representation of WA n -i, i.e. by 
(6.2.2) with C no -i = C n -i = A n -i and p = n + 1, q = n + 2. In particular, the central 

charge for a Z n -model equals c = 2 "+2^ • ■'■ n ^ ms sec ti° n we have explicitly verified for 
n = 3 and 4 that the representations of the twisted sector of WA n -i correspond to the 
spectrum of the Z n -model with twisted boundary conditions. In fact, this is also true for 
the Zs-version [65]. These explicit results are in agreement with the statement that the 
field content of the spin quantum chain (3.1.1) at (j) = ip = 0, A = 1 with twisted boundary 
conditions Tn + i = can be described by a representation of a twisted WA n -i for all n. 
Thus, it is possible to calculate the spectrum of the twisted Z n -quantum chain by (6.2.2) 
using C n -i = An-i, £ no = C[nj and p = n+l, q = n + 2. 

For cyclic boundary conditions (rjv + i = cu _fi Ti, < R < n) the diagonal symmetry of 
the statistical mechanics model is known to be broken such that the spin h — h takes on 
rational values. The dimensions of the chiral fields, however, are unaffected by this change 
of boundary conditions. This has been verified in [37] and [71] for the case n = 3, in [14] for 
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n = 4 and more abstractly for general n in [77]. In [37] a similar result has been obtained 
for r^v+i = u~ R r^ , 0<_R<3, n = 3 and the only effect of a factor u~ R for all n should 
be to combine the left- and right-chiral parts in a non-diagonal way. 

We would like to conclude by mentioning that the case of free boundary conditions should 
be described by the representations of the orbifold of WA n -i at its first unitary minimal 
model. This can be read off from the one-to-one correspondence of boundary conditions of 
W-algebras and spin chains which we have just observed. After mapping the cylinder on 
which the CFT lives to a strip in the plane with open ends, only those fields can survive 
which are invariant with respect to transformations at the boundary which implies that 
precisely the orbifold is expected to survive in the case of free boundary conditions. 

6.4. Classical W-algebras and Virasoro structure 

In this section we first give a brief review of the notion of classical W-algebras - for details 
see e.g. [59,60]. Afterwards, we make some general statements following from Virasoro 
covariance which are only partially known in the literature. 

For a characterizations of a classical W-algebra we closely follow the quantum case - 
compare section 6.1: 

1') A classical W algebra is a vector space of local chiral classical fields, i.e. there is a 
basis of fields (f>(z) depending on a complex coordinate zeC with definite ' conformal 
dimension' d(<fi) and mode expansion 

cf>(z) = ^"^Vn (6.4.1) 

n-d(»eZ 

such that all fields <p(z) are either bosonic or fermionic and that their dimensions are 
half-integer 2d(</>) E Z. 
2') These fields form a Poisson bracket algebra which we write in OPE form. The 'classical 
OPE 1 of two fields 4>{z), x( w ) is °f the form: 

d(<P) + d( x )>d(i,) 

The are the structure constants of the classical W-algebra. Instead of the Wick 
rules of the the quantum case, the classical OPE has the derivation property of Poisson 
brackets. For example, {4>(z), xi™)^^)} = V'C 10 )} + {4>{ z )-> xiw)}^^)- 

Inserting the mode expansion (6.4.1) into (6.4.2) and performing Cauchy integrals, 
the classical OPE gives rise to a Lie bracket - like on the quantum level. 
3') The classical W-algebra contains a distinguished field, the energy-momentum tensor 
L with classical OPE: 

rw \ ti m c 2L(w) L'(w) .„ A , 

{L(z), L(w)} = ^ + ^ + -^-^ + reg. 6.4.3 

2 (z — wr [z — wY (z — w) 
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4') 



The energy-momentum tensor L defines a grading on the classical W-algebra: 



d(<l>)<f>(w) ^ At^(w) 

d(i,)^d(4,) 
2+d(4>)>d(4,) 

where d((p) is the conformal dimension of the classical field <fi. 
5') Classical fields transforming as 

/rn u dj^jw) : <j>'{w) ^ A^ipjw) 

{L(z),(p(w)} = -. -7T + - + > M + reg. (6.4.5) 

W V /J (z-«;) 2 (z-w) ^ (z - ty)2+dW-d(V) y v y 

d(<A)>ci(V)+i 

are called 'quasi-primary'. The field is called 'primary' if the remaining sum is 
absent. 

6') All classical fields in the classical W-algebra are differential polynomials (with finite 
order) in a countable basic set of quasi-primary fields {W&} - the so-called ' gener- 
ators' 1 In general, this differential polynomial ring will not be freely generated 
but may have relations, i.e. differential polynomials that vanish identically. Note that 
in contrast to the normal ordered product of the quantum level, the product of two 
classical fields is (anti-)commutative. 

7') The 'phase space" 1 of the classical W-algebra is the vector space of meromorphic func- 
tions from the Riemannian sphere C into the complex numbers C. The classical fields 
are defined as integral kernels of functionals on this phase space. In particular, all 
classical OPEs (6.4.2) have to be understood after Cauchy integration with some 
meromorphic test function. 

Note that usually a notation different from the above is used in the literature. Usually one 
uses the compact manifold S 1 instead of the Riemannian sphere C. Then the poles in the 
OPE are replaced by 5-functions on the r.h.s. of the Poisson brackets. Test functions are 
smooth square-integrable functions on S , not meromorphic functions. However, this is 
only a minor difference in the point of view since after all we are interested in the algebraic 
structure of a classical W-algebra, and here both formulations lead to the same results. 

The description in section 6.1 of Kac-Moody algebras, cosets, orbifolds and DS reduction 
apply at the classical level as well with some obvious minor modifications. Therefore we 
do not explicitly repeat the corresponding discussion of section 6.1. 

In contrast to the quantum case one does not have a bilinear form on a classical W- 
algebra which one can use to prove basic theorem e.g. about the su(l, 1) structure of the 
space of fields. Nevertheless, explicit constructions are known on the quantum level (see 
e.g. [30]) which can be transferred to classical W-algebras. Then one can use these explicit 
constructions in order to prove e.g. that the generators of a classical W-algebra can be 
chosen primary in quite general situations. In the remainder of this section we wish to 
indicate how this should be done. 

A first result is a formula for a quasi-primary projection of two classical fields including 
some derivatives: 
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Lemma II: ( Classical quasi-primary product) For any two quasi-primary fields 4>, x let 
e<"(«,x):= £ (-l)-( M( * ) J; r " 1 )( MW ^ r " 1 )(8**)(8"x). (6A6) 

is a quasi-primary field of conformal dimension d(<f>) + d(x) + v. ■ 

In order to prove this statement one has to apply (6.4.4) to the r.h.s. of (6.4.6) and check 
that it has the correct transformation properties. Of course, an equivalent expression can 
be obtained by applying a suitable classical limit [28] to the formula for the quantum case 
[30]. Covariance problems of this type have been treated in the mathematical literature 
already many years ago. For example, finding the projection of the rth derivative of 
the product of two modular forms onto a modular form (see e.g. [128]) is equivalent to 
determining the quasi-primary projection. In fact, we have taken Eq. (6.4.6) from the 
formula (1) of [128] which is called the 'rth Rankin-Cohen bracket'. 

Note that it is possible that the classical quasi-primary product (6.4.6) vanishes identically. 

Now we turn to the question of Virasoro covariance of the space of fields of a classical 
W-algebra. To this end we need a weak property of the generating set. Then we can show 
by construction that for any W-algebra with this property the space of fields decomposes 
into quasi-primary fields and derivatives thereof. We will also say that a W-algebra has a 
quasi-primary basis if it admits such a decomposition. 

Def. VII: A classical W-algebra is called 'regularly generated' if for each fixed dimension A 
there are only finitely many generators with conformal dimension less than A. 

Claim I: Any classical W-algebra that is regularly generated admits a basis of quasi- 
primary fields and derivatives thereof. ■ 

This statement should be proven for homogeneous polynomials of a fixed degree in the 
generators. For the linear polynomials the claim is true by definition. Also for the second 
order polynomials it is immediately clear that the formula (6.4.6) yields precisely as many 
linearly independent quasi-primary fields as are needed in order to complement the deriva- 
tives to a basis. More precisely, all linear relations among the second order quasi-primary 
fields are given by Q^ r \(j),x) = (%, (f>) and Q^ r '{4>-i 4>) = for bosonic and r odd 

or <fi fermionic and r even. For polynomials of higher degree we forget at the moment 
about the relations possibly present in the space of fields (they can be imposed at the 
very end). Then one has to gain control about the action of the derivative. This can e.g. 
be done by using partitions into different colours (compare [62] for a presentation in the 
quantum case). Now a basis of the quasi- primary fields can be labeled by these partitions. 
It remains to be checked that those quasi-primary polynomials associated to the partitions 
are indeed linearly independent. We leave this inspection for future investigations. 

Note that this claim should in particular be valid for W-algebras that are generated by 
finitely many fields whose dimensions are strictly positive. However, the above argument 
should apply to regularly generated classical W-algebras without additional effort. Claim 
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I might be true in even more general situations but this will certainly give rise to extra 
complications in the proof. 



If a classical W-algebra admits a basis of quasi-primary fields and derivatives thereof one 
can derive a specialization of the classical OPE (6.4.2) when applied to quasi-primary 
fields: 

Lemma III: ([61]) For any classical W-algebra with a quasi-primary basis the Poisson 
brackets of two classical quasi-primary fields (f>(z), x( w ) have the form 



{<t>(z), x (w)}= _ w)dW+d(x) + E c; x ^ 

^ ' ip quasi-primary 'p — Q 

d(<#>) + d(x)>d(V) (6.4.7) 

(d(<t>)-d(x)+d(*l>)+r-l^Qr^( w j 

r!( 2dW+r_1 )(2 - w yW+d{ x )-d{ip)-r + reg - 

m 

This statement is also true at the quantum level (see e.g. [30]), and in fact was first 
discovered there. It can be proven along different lines. One can either use the realization 
of L_i as the differential operator L_\(j>(z) = z(zd + 2d(<p))(/)(z) acting on quasi-primary 
fields <j>(z). Then one uses covariance of (6.4.2) with respect to this differential operator 
[61]. Equivalently, lemma III can be derived from the Jacobi identity for the Poisson 
brackets {L(y) : {</>(z), x( w )}} + cyclic = using (6.4.5). 

In view of lemma III we can introduce the following shorthand notation for the Poisson 
brackets of two classical quasi-primary fields 0, x writing only the quasi-primary fields on 
the r.h.s.: 

ip quasi-primary 
dO) + d(x)>d(V0 

The notation (6.4.8) has to be understood in the sense that the arguments and poles 
have to be re-instituted in the form of (6.4.7). Similarly, the coupling constants involving 
derivatives can be recovered from (6.4.7) since they are uniquely given by those involving 
only quasi-primary fields. In this shorthand notation the classical Virasoro algebra (6.4.3) 
becomes L* L = | + 2L. 

This shorthand notation is useful for arguing that one can choose the generators of a 
classical W-algebra primary in a quite general situation. 

Claim II: All generators of a classical W-algebra except for the Virasoro field L can 
be chosen to be primary if the W-algebra is regularly generated, the Virasoro centre is 
non-zero c ^ and the conformal dimensions of the generators are strictly positive. ■ 

A sketch of the proof of this statement is the following: If claim I is true there is a 
quasi-primary basis according to the assumptions and we can use the shorthand notation 
(6.4.8). First we note that the generators W^> of conformal dimension d(W^') < 2 
must already be primary because of the assumptions. Now consider any generator W^> of 



d(<f>)+d(x)-d(il>)- 
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conformal dimension d(W^) = 2. Since W^> is quasi-primary, the OPE with L must read 
L * WW = ajc + 2W {j \ Then = - 2a jL is primary. For all other dimensions 
we first observe that 

L * flM (L , = + r - ffc+i>!, + W) + r + 2) eM (L, ,) 

+ J] «},*fi , ' ) (i,l» 

?/> , s 

d(i/i) + s + 2<d(c#>)+r 

The ^ are certain structure constants. For the generators of conformal dimension 
not equal to two, the OPE with L must have the form 

L*W® = d(W (j) )W^ + a ^ (fc) (6.4.10) 

k<d(WU))-2 

with one quasi-primary field per conformal dimension k. Now we can use (6.4.9) to 
add fields (L, (p^) with suitable parameters. Since there may be further fields ip of 
the same dimensions on the r.h.s. of (6.4.9) this procedure may have to be iterated adding 
further fields with suitable parameters. According to the assumptions there are only finitely 
many fields below a fixed conformal dimension. Therefore, after adding finitely many 
counterterms one is left with a linear system of equations that must be solved. Once this 
system is solved, the generator becomes primary and claim II is verified. From (6.4.9) one 
obviously needs c 7^ in order to have an invertible system. For low conformal dimensions 
(e.g. d(W^) < 4) one can argue that this system is indeed invertible. The general case 
may become quite complicated and we will therefore leave the proof of invertibility to 
future investigations. 

We would like to conclude this section by remarking that from here on it is in principle 
straightforward to construct classical W-algebras by making some assumptions e.g. on the 
conformal dimensions and then checking Jacobi identities. For example, one can check 
that there is a unique classical W(2, 3) (up to isomorphism) if the central terms are chosen 
to be non-zero. 

6.5. Classical cosets and orbifolds 



In this section we discuss orbifolds and cosets of classical W-algebras. During the attempts 
to understand a mysterious algebra of type W(2,4, 6) [51] using some classical limit evi- 
dence emerged that a classical analogon would need more generators than this quantum 
W-algebra and that the same statement should be true for the orbifold of of the N = 1 
quantum super Virasoro algebra (for more details see [28]). It was then argued in [31] 
that generic classical cosets and orbifolds are -in contrast to the quantum level- infinitely 
generated and satisfy infinitely many constraints. Let us now briefly summarize the results 
of [31]. 

Classical orbifolds are defined in the same manner as they are defined at the quantum 
level (see the end of section 6.1). Below we restrict to Z 2 automorphisms p that act on the 

63 



(6.4.9) 



finitely many generators {W a | a G X U /C} as follows: 

p(W a )=W a VaG/C, 

(6.5.1) 

p (W h ) = -W b Wbel . 

We can divide the index set X into two subsets: A set Xi referring to bosonic fields and 
a set X2 referring to fermionic fields transforming nontrivially under the automorphism p. 
It is easy to determine a generating set classically. Note that the nontrivial p-invariant 
differential polynomials are even order in the {Wf, \ b G X}. Plainly, every even order poly- 
nomial can be regarded as a polynomial in quadratic expressions. Therefore the quadratic 
expressions formed out of the {W b | b G X} generate the orbifold together with the invariant 
fields {W a I a G /C}. A redundant set of quadratic generators is given by: 

X l b >i := W b (t) W^ b,cel, < !,j £ Z (6.5.2) 

where := d l W b - The derivative acts on the generators (6.5.2) as follows: 

dXi : i = X^ + Xif\ (6.5.3) 

Using the action of the derivative (6.5.3) and paying attention to the Pauli principle for 
the fermionic generators, i.e. that fermions have odd Grassmann parity, one can choose 
the following minimal set of generators for the orbifold: 



W a , a G /C (invariant generators) , 



X° h i := W b d J W c , b < c, 6, c G X, < j G Z , 



a, a 



= Wdd^Wd, d G Xl, < j G 2Z (square of bosons) , 

= W e d^W ei e G X 2 , < j G 2Z + 1 (square of fermions) 



(6.5.4) 



where '6 < c' denotes some ordering of the original generators. Eq. (6.5.4) shows that Z2 
orbifolds are always infinitely generated at the classical level (if there is one field trans- 
forming non-trivially under p). 

In order to find the complete set of relations we first regard all W a as independent. The 
complete set of relations satisfied by the redundant set of generators (6.5.2) is generated 
by 

K% ~ e b ,cXH = 0, (6.5.5a) 
Xli X di - ^ d K' k d X J d = 0, (6.5.5b) 

where e ^ c = — 1 if both W b and W c are fermions, and e ^ c = 1 otherwise. (Clearly, choosing 
certain indices in (6.5.5) equal leads to trivial relations). The proof that (6.5.5) indeed gen- 
erate all relations is a simple sorting argument. Plainly, a basis for invariant polynomials is 
given by those differential monomials that have even degree in the generators transforming 
non-trivially under p if the generators are arranged in a suitable order. Obviously, any 
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given monomial in the fields (6.5.2) can be re-sorted into this order using the relations 
(6.5.5) and (anti-)commutativity of any two fields. 

It is straightforward to derive the relations satisfied by the nonredundant set of generators 
(6.5.4) from (6.5.5). One simply has to recursively apply (6.5.3) (which encodes the action 
of the derivative) in order to express the relations (6.5.5) in terms of the generators (6.5.4). 

Next, we further elaborate some of these relations for two examples. One of the 
simplest examples of orbifolds is the bosonic projection of the N = 1 super Virasoro 
algebra. The N = 1 super Virasoro algebra is the extension of the Virasoro algebra L by 
a primary dimension | fermion G. According to (6.5.4) a nonredundant set of generators 
for the classical orbifold is 

L, $ n := Gd n G for all odd n. (6.5.6) 

In particular, this orbifold has one generator at each positive even scale dimension. Using 
the notation of (6.5.4) we have the identification $ n = X 0,n if we omit the irrelevant lower 
indices. From (6.5.5b) one reads off = X°>iX°> 1 +X°'°X^ 1 = X°^X ' 1 because X ' = 0. 
In terms of the generators (6.5.6) these infinitely many relations read 

$ n $ m = q for any o < n, m G 2Z + 1. (6.5.7) 

In this case the particular subset (6.5.7) of relations for the nonredundant set of generators 
can also immediately be inferred from the Pauli-principle: $ n $ m = (Gd n G)(Gd m G) = 
—G 2 (d n G)(d m G) = 0. However, (6.5.5) encodes more relations. For example (6.5.5) and 
(6.5.3) imply X ' 1 ^ ' 1 = X Q '\X 2 ^ + 2X 1 ' 2 + X ' 3 ) = X^X 1 ' 2 = -X^X 1 ' 2 . For the 
nonredundant set of generators this implies the following relation at scale dimension 10: 

^d 2 ® 1 = 0. (6.5.8) 

As a simple second example we consider two commuting copies of the Virasoro algebra 
(Li and L 2 ) with equal central charges. Then W := L\ — L 2 is primary with respect to 
L := Li + L 2 . Furthermore, p(L) = L and p{W) = —W is an automorphism of this 
W(2, 2). According to (6.5.4) the subspace invariant under p is generated by the following 
fields: 

L, $ n := Wd n W for all even n. (6.5.9) 

Again, we obtain one generator at each positive even scale dimension. In this case rewriting 
the relations (6.5.5) for the redundant set (6.5.3) in terms of the nonredundant set (6.5.4) 
is slightly more complicated. Using the notation X 8J = W'''W^ one checks that from 
(6.5.3) and (6.5.5a) d 2 (X°> X°> ) = SX ' 1 ^ ' 1 + 2X >°d 2 X°> and d 2 X°>° = 2X 1 ' 1 + 
2X°> 2 . Using (6.5.5b) it is straightforward to check that d 2 (X°> X°> ) + 6X >°d 2 X°> + 
8X 0,0 X 0,2 = 0. This relation arises at scale dimension 10, which is the lowest scale 
dimension admitting a relation. In terms of the generators (6.5.9) it reads 

a 2 ($°$°) - 6$°a 2 $° + 8$°$ 2 = 0. (6.5.10) 

Now we turn to a discussion of cosets. First, we briefly discuss cosets of the type 
W/Q with Q as simple Lie algebra. By definition, the generating set of the coset is a 
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generating set for those differential polynomials in W that are invariant with respect to 
the action of Q. Problems of this type are known in the mathematical literature under the 
name of ' invariant theory' and one can use their results (see e.g. [124]) in order to obtain a 
generating set. Observing that Qk C W and that there is always an invariant bilinear form 
on Q it is immediately clear that the coset will have infinitely many generators because 
any bilinear expression in the currents of Qk with an arbitrary number of derivatives is an 
invariant field. Concerning the relations one also sees that there will be infinitely many 
of them because non-trivial invariant polynomials in finitely many variables always satisfy 
some relations. Now one can put infinitely many derivatives (like for the generators) and 
it is clear that there will be infinitely many relations. 

Cosets of the type W/Qk are a bit more complicated. The following result was obtained in 
[31]: Assume that the restriction of the central term of W to Qk is nondegenerate (we do 
note assume Q to be simple) and that it is possible to partition the generating fields of W 
into the generating fields J a of Qk and a complementary set of generating fields that 
form primary multiplets with respect to the current algebra Qk (to avoid confusion note 
that the Jf- need not have conformal dimension one). Under these conditions the adjoint 
action of Qk on its complement W is a representation of Q: 

{J a (z), JHw)} = £ V w ( 6 - 5 - n ) 

3 

in terms of matrices T(a). Now one introduces a covariant derivative V acting on the fields 

VJt := dJt -\Y, Jb9 ha T(a)lJj- (6.5.12) 

a,b,j 

where 

g ba Jg 

the inverse of the metric on Q. Note that the covariant derivative T> is indeed 
a derivation. The reason simply is that if a field <j> transforms under a representation T 
and x under a representation T then 4>d n x transforms under the tensor product of the 
representations T and T. 

Now the coset W/ Qk is the space of those ^-invariant differential polynomials in the fields 
j/- where as derivative the covariant derivative (6.5.12) has to be used. From this charac- 
terization it is immediately clear that Wj Qk must have infinitely many generators satisfying 
infinitely many constraints precisely for the same reasons as for W/Q. 

A simple consequence of this general result is that for classical cosets of type Qk ® 
'Hm/'Hik+m (where i denotes the Dynkin index of the embedding Ti C Q) one has 

lim Qk®n m = Qk (6 513) 

— 7i lk+ m h 

and that the coset on the l.h.s. for generic m is a deformation of the coset on the r.h.s. 
Eq. (6.5.13) is also true for cosets of quantum W-algebras (see e.g. [34]). 



66 



Let us now conclude this general discussion of classical cosets and turn to some more 
concrete statements. First we observe that if one wants to perform quasi-primary pro- 
jections in coset W-algebras one may use the formula (6.4.6) if one replaces the ordinary 
derivative d by the covariant derivative T>. The reason is mainly that T> — d commutes 
with the coset energy-momentum tensor L. 

For cosets of Kac-Moody algebras the coset energy-momentum tensor L is always com- 
posite which implies c = 0. Then L is obviously primary. Under these conditions the 
following fields are primary as well if 4> and x are primary fields of dimension d(4>) and 
d(x) respectively: 

Q (0) (</>,%) = <I>X, ±Q {1) (<P,x)=d(<j>W -d( X )<P'x- (6-5.14) 

As a curiosity aside we note that the covariance properties of the second expression have 
already been noted several decades ago in a different context (see e.g. eq. (24) on page 23 
of [125]). It should be clear to the reader that the assumptions of claim II are not satisfied 
for composite L and thus there need not be a primary generating set. 

Let us now illustrate these basic remarks with a few statements on some classical cosets. 
The absence of a primary generating set can be seen very nicely in the classical coset 
sl(2)/U(l) (compare also [122]). The fields 

W n+2 := Q {n) (J + , J~), 0<neZ (6.5.15) 

form a quasi-primary generating set of the classical coset sl(2)/U(l) (for the relations see 
[31]). The fields (6.5.15) cannot have a primary projection unless they are already primary. 
It can be checked explicitly that Wi = kL and W3 are primary whereas W4, W$ and Wq 
are not primary. This demonstrates that there is no primary basis in the classical coset 

sU2)/mi) 7). 

Also the classical coset SVIR(N = 2)/U(l) can be controlled quite well. One finds pre- 
cisely one generator for each integer scale dimension greater or equal to two. For this 
coset the energy-momentum tensor L has a non-vanishing central term because it is not 
composite. Thus, according to claim II it should be possible to choose the generators with 
dimension greater than two such that they are primary. This can in fact be checked explic- 
itly for small conformal dimensions (we have performed checks up to conformal dimension 
five). For more details on this coset see [31]. 

Finally, let us briefly consider the classical coset sl(2) k / sl(2). From the discussion of 

the (3 — 7-system in [31] one expects that the classical coset sl(2) k /sl(2) should have a 
specific Poisson subalgebra at k = — \. In order to check this one can perform explicit 

computations of Poisson brackets in sl(2) k /sl(2) at low scale dimension. We use normal- 
ization conventions analogous to [28]. According to [31] we have to consider the following 

7 ) Note, however, that due to (6.5.14) there are more primary fields in this coset than 
claimed in [122]. 
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generators: 



S m , n : = (d m J + )d n J- + ±(<9 m J°)d n J° + (d m J-)d n J+ , 
S m ,n,k ■■ = (d m J-)(d n J°)d k J + - (d m J-)(d n J + )d k J° - (d m J°)(d n J-)d k J + (6.5.16) 
+ (d m J°)(d n J+)d k J- + (d m J + )(d n J-)d k J° - (d m J + )(d n J°)d k J- . 

Eliminating the redundant generators from (6.5.16) we see that this classical coset is of type 
W(2, 4, 6, 6, 8, 8, 9, 10, 10, . . .). We can choose the following non-redundant quasi-primary 
generators up to scale dimension 6: 

L = 2k S 0fi , 

W4 = S 0; 2 + 2^1,1 i 
^ = 50,4-1051,3+1052,2, 

Web = 5o,i,2- 

One can check that L and Wq^ are primary whereas W4 and Wq & are only quasi-primary 
which once again demonstrates the absence of a primary generating set. Looking for a 
Poisson subring one finds the following useful linear combination of generators: 

W 6 := ^W 6a -2W 6h . (6.5.18) 

We have checked up to scale dimension 6 that L, W4 and Wq close among themselves under 
Poisson brackets for all k (i.e. we have not looked at any fields beyond scale dimension 6 
on the r.h.s. of the classical OPEs - still a quite non-trivial fact). From this observation 
and the study of the (3 — 7-system in [31] it seems sensible to expect that the classical coset 
sl(2) k /sl(2) has a Poisson subalgebra generated by specific even dimensional fields for all 
values of the level k. 



(6.5.17) 



Concerning the quantum case of orbifolds and cosets in general there is the following 
hypothesis: Upon normal ordering of the classical relations correction terms emerge that 
contain classical generators such that one can eliminate precisely one generator per genera- 
tor of the classical relations. This 'cancellation' has been checked so far rigorously only for 
a very simple coset [31] and for a few dimensions in some other cosets and orbifolds [31, 28]. 
Note that the character argument which has been used in the argumentation supporting 
the results of the next section relies heavily on the assumption that this cancellation always 
happens. 

6.6. The structure of quantum W-algebras 

One of the main aims of the explicit constructions of quantum W-algebras e.g. in refs. 
[30, 102, 94, 51] was to gain some insight into the structure of quantum W-algebras and 
the associated rational conformal field theories (RCFTs). The rational models of those W- 
algebras existing for isolated c only can be either interpreted as extensions or truncations of 
W-algebras obtained by DS reduction, or the effective central charge is integer (see [49] for 
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a summary) . After identifying the Casimir W-algebras and the orbifold of the N = 1 super 
Virasoro algebra among the W-algebras existing for generic c two unexplained solutions 
remained: An algebra of type W(2,4, 6) [102] and an algebra of type W(2,3,4, 5) [94] 
which both give rise to generic null fields. In particular the origin of this W(2,4, 6) was 
unexplained for some time [51]. In the meantime both of these algebras have been identified 
[31, 27, 28] and it was realized that they fit in a general new structure: They play the role 
of 'unifying W-algebras'' for Casimir W-algebras [27,28]. The purpose of this section is to 
explain this new unifying structure. 

One of the basic observations is that truncations of Casimir W-algebras are quite general, 
i.e. that for certain values of the central charge c some of the generators turn out to be 
null fields. This can be investigated systematically using the Kac determinant [33] for the 
vacuum representation. This has been carried out in [27] leading to the truncations in 
table 8. 



Casimir algebra 


c 


truncated algebra 




CA n (n + r,n + s) 


W(2,...,rs- 1) 




cs„(2n-2 + r, 2n - 1 + s) 


W(2, 4, . . . , rs — 2) for r-s even 


wc n 


c Cn (n + r, 2n - 1 + s) 


W(2, 4, . . . , rs — 2) for r-s even 


Orb(WP n ) 


c-v n (2n - 3 + r, 2n - 3 + s) 


W(2, 4, . . . , rs — 2) for r-s even 



Table 8: Truncations of Casimir W-algebras 



The parameterization of the central charge in table 8 is given by the following generalization 
of (6.2.2): 

c c J P , q) =n-n { MZlpL (6.6.!) 

where p v is the dual Weyl vector (for simply laced C n one has p v = p). 

The result of table 8 gives rise to the structure indicated for the Casimir W-algebras WA n 
in Fig. 4. The horizontal lines correspond to the algebras WA n that exist for generic 
central charge c (or equivalently, generic 'level' k). For the kth unitary minimal model all 
generators of WA n with dimension greater than k 2 + 3k + 1 become null fields leading to a 
truncation of WA n to an algebra of type W(2, . . . , k 2 + 3k + 1) for all sufficiently large n. 
It was argued in [27] that one can interpolate these truncations at fixed k to non-integer 
n giving rise to a new W-algebra existing for generic c - a 'unifying W '-algebra' . 

In general we will call a W-algebra a unifying W-algebra if it is finitely generated, exists 
for generic c and has the property that a discrete family of W-algebras exists (e.g. Casimir 
W-algebras based on a Lie algebra of fixed type) such that infinitely many members of 
this family at certain distinct values of the central charge truncate to this unifying algebra. 
Since unifying W-algebras exist for generic c they can e.g. be thought of as continuations 
of Casimir W-algebras W£ n to real values of the rank n for certain values of the cen- 
tral charge. Usually, unifying W-algebras give rise to null fields, i.e. they are non-freely 
generated. 
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The part of Fig. 4 corresponding to unitary minimal models (k a positive integer) can be 
understood in terms of level-rank duality [9, 34] using the coset realization of WAn-i' 



WA n -! = 



sl(n) k esl(n) 1 _ sl(k+l) 



sl(n) 



k+i 



sl(k) n ®U(l) 



= CV{k) . 



(6.6.2) 



n 

9 
8 
7 
6 
5 
4 
3 
2 
1 



W(2,3,4,5) 
U(l) 



'10 



Zc 



1 



H h 

2 3 



fc e 



6 



(c generic) 
I 

W(2, . . . , k 2 + 3k + 1) 
(at level k) 



WA 2 = W(2,3) 
ytMi Vir 



7 



8 9 



k 



Fig. 4: The structure of the Casimir W-algebras WA n and their unifying 
W-algebras. 



More precisely, it was shown [9, 34] that the cosets in (6.6.2) give rise to equivalent energy- 
momentum tensors. In [28] evidence was collected that (6.6.2) is in fact true on the level 
of W-algebras. Thus, the symmetry algebra of the CV(k) model is a unifying W-algebra 
for the kth unitary minimal model of WA n -i- Note that the l.h.s. of (6.6.2) is defined for 
integer n and arbitrary k, whereas the r.h.s. is defined for integer k and general n. The 
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isomorphism in (6.6.2) is valid iff k and n are both positive integers. Note that CV(k) 
does indeed exist for generic c or fc, thus giving rise to a coset realization of the unifying 
W-algebra predicted by inspection of the Kac-determinant. 

It was argued in [28] that the field content of CV(k) is given by 

CV(k) = W(2, 3, . . . , k 2 + 3k + 1) (6.6.3) 

with two generic null fields at dimension k 2 + 3k + 4 which is consistent with the prediction 
from the Kac-determinant. 

Now we would like to discuss the first unitary minimal models of the algebras WA n -i in 
some more detail. On the one hand, the case k = 1 corresponds to the Z n parafermions 
which are relevant to the Z n -quantum spin chains that we have discussed earlier. On the 
other hand, the coset sl(2) n /U(l) arising at k = 1 obviously has two effective degrees of 
freedom, and using the character argument of P. Bouwknegt (see e.g. [33]) it is obvious that 
it should be finitely generated. Nevertheless, there is some confusion about this question 
of the generating set in the literature (see e.g. [18, 112, 127]). Therefore we would like to 
explicitly recall what we have been able to show with some effort in [28]: 

Proposition III: ([28]) The quantum coset sl(2) n /U(l) has a closed W(2, 3,4, 5) subal- 
gebra with two generic null fields at conformal dimension 8. Furthermore, there are no 
additional generators of conformal dimension 6, 7 or 8 in the coset sl{2) n /U{l). ■ 

Remarks: 

1 ) More strongly, a character argument [28] indicates that the coset sl(2) n /U(l) is iso- 
morphic to the W(2,3,4,5). Unfortunately, we have not been able to prove this 
stronger statement. 

2) This W(2, 3,4, 5) is the previously unexplained algebra found in [94]. 

We have already mentioned that this W(2, 3, 4, 5) unifies the first unitary minimal models of 
WAn-i that are related to critical points in the Z n -spin quantum chains of earlier sections. 
In view of sections 6.2 and 6.3 we would be interested in the twisted sector of this algebra. 
Indeed, an inner automorphism of sl(2) n gives rise to an outer Z2 automorphism of the 

coset sl(2) n /U(l) [28]. Unfortunately, at the moment we do not know how to apply the 
standard machinery of coset representations in the presence of an inner automorphism 
which means that in contrast to the periodic sector the twisted sector cannot be easily 
obtained from the coset construction. Nevertheless, one can argue that the orbifold is of 
type W(2, 4, 6, 8, 10) with a generic null field at scale dimension 14 for generic n. According 
to section 6.3 this orbifold should be relevant for the Z n -spin chain with free boundary 
conditions. We conclude this discussion by recalling the simple observation that due to 
the coset realization the irreducible modules of Z n parafermionic model can grow at most 
as fast as a free module of two chiral bosonic fields. It would be interesting to see if the 
presence of a unifying object at generic real n can be used for obtaining more insights into 
Z n parafermionic conformal field theories. 

The investigation of the Kac determinant is not restricted to integer k but also applies 
to k G Q - including in particular the non-unitary minimal models. Using the parameteri- 
zation (6.6.1) the field content can be read off from table 8. As before, one can interpolate 
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these further unifying W-algebras to generic n (or c) . These unifying W-algebras can also 
be realized in terms of cosets. Now, one needs Drinfeld-Sokolov reductions based on sl(r) 

for a non-principal embedding such that they have a sl(k) © U(l) Kac-Moody subalgebra: 

yvf% k = W(l fc2 , 2, 3, . . . ,r - k, (^) 2k ). (6.6.4) 

The k 2 currents form a sl(k) © £7(1) Kac-Moody algebra, the fields of dimension 2, . . ., 
r — k are singlets with respect to this Kac-Moody and the 2k fields of dimension r ~^ +1 
are a Kac-Moody multiplet. More precisely, these Ik fields transform as two £7(l)-charge 
conjugate defining representations of sl(k). In [28] it was argued that 

WAn-i = r k, l^ = W(2,3,...,(/c+l)r + fc) at c A _ 1 (n + /c,n + r) (6.6.5) 

sl(k)®U(l) 

with two generic null fields at dimension (k + l)r + k + 3. Identifying W^+i +1 ^ with the 

unconstrained sl(r + 1) Kac-Moody algebra, the case r = k + 1 is identical to the CV(k) 
cosets which we have already discussed. One can also apply (6.6.5) to the case k = 0. In 
this case, the sl(2) embedding is the principal one and there is no Kac-Moody subalgebra. 
Thus, the coset in (6.6.5) is just WA r -i and we recover eq. (2.4) of ref. [27] (note that 
this equality was originally observed in [106]). 

The Casimir algebras WA n in Fig. 4 obviously exist also for negative k and from the 
preceding discussion we conclude that one can make sense of the unifying W-algebras also 
for negative n. However, it is also possible to directly make sense of WA n at negative n 
[96] (compare also the comments on WD- m below) or even for complex n G C [53]. Among 
those, the WA n at negative n are defined less formally since they can be given a coset 
realization WA- n -i = (sl(n) k © sl(n) _ t ) / sl(n) k _ 1 and have a definite finite generating 
set [96]. 

Using WA n for real n or the fact that Fig. 4 is densely covered by unifying W-algebras one 
can define [96] a l universaF W-algebra for all Casimir algebra based on A n that depends 
on two parameters, e.g. k and n. For generic irrational k ^ Q this universal W-algebra will 
have infinitely many generators. Doing so, one can nicely fit the Woo and Wi+oo algebras 
(see e.g. [17]) in the picture of Fig. 4: They can be recovered from the topmost border 
n = oo respective right border k = oo. Then, e.g. the truncations of Wi+oo at c = n [64] 
arise naturally from Fig. 4 (for more details see [28]). 

To summarize, the picture of the space of all Casimir algebras WA n and their unifying 
W-algebras as presented in Fig. 4 is mainly a conjecture. However, this conjecture has 
been thoroughly checked in [27,28]: Beyond the inspection of the Kac determinant and 
cross-checks with limiting cases which we mentioned above one can perform character 
computations for the cosets and compare structure constants or minimal models (where 
known). Even more, the simplest of the relevant cosets can be constructed explicitly [28] 
and e.g. the statement about the generating set can be quite rigorously verified. 

The structure of the Casimir algebras WV n: WC n and W£> n is similar to Fig. 4 which 
we have discussed in some detail. The truncations inferred from the Kac determinant are 
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presented in table 8. Whereas in the case of WA n we were able to provide a coset realization 
for all unifying W-algebras, coset realizations of unifying W-algebras for the other Casimir 
algebras are only partially known. The known coset realizations for unifying W-algebras of 
Casimir algebras are summarized in table 9 (for a detailed discussion see [28]). Sometimes 
an additional Z 2 orbifold is needed which is denoted by 'Orb' in table 9. If 'Orb' is put in 
brackets only part of the algebras give rise to non-trivial outer automorphisms - for the 
other ones we do not have to take an orbifold. 



Casimir 
algebra 


central charge 
c 


coset realization 
of unifying algebra 


dimensions of 
simple fields 


dimension of 
first null field 




c An ^ 1 (n+k,n+r) 


s 


r-k,l k 
l(k)(£U(l) 


2, 3, . . . , kr + r + k 


kr + r + k + 3 




c Bn (2n+k-l,2n+l) 
c Bn (2n, 2n+k) 


(Orb) ( 


/ so(fc) K © S o(fc) 1 \ 


2,4, ...,2k 


2A; + 4 


wc n 


cc n (n+k+l,2n+2k+l) 


S p(2fc)„_i 


2, 4, . . . , 2k 2 + 4k 


2k 2 + 4k + 5 


Ovb(WV n ) 


c Vn (2n+k-2, 2n+k-l) 


Orb 


f S o(fc+l) 2 \ 

V s °w 2n J 


2,4, ...,/c 2 + 3/c 


k 2 + 3k + 4 



Table 9: Coset realizations of unifying W-algebras 



Among the unifying W-algebras for these three families of Casimir algebras those for WC n 
are particularly interesting because they can in some sense be regarded as algebras of type 
'WD-fc' - continuations of the WV m series of Casimir algebras to negative m. Historically, 
the mysterious W(2,4, 6) was the starting point for the discovery of unifying W-algebras 
in [27,28]. It was proposed in [95] that this W(2,4, 6) can in some sense be regarded as 
belonging to the series of WE> n Casimir algebras: W(2,4,6) = WP-i. Furthermore, all 
known minimal models of this W(2, 4, 6) [51] were identical with minimal models of some 
WC n Casimir algebra. Only after generalizing these observations we became aware that 
the algebras WV_ k can be given a group theoretical meaning. 

Negative dimensional groups are known in the literature (see e.g. [46]). More precisely, 
some aspects of the representation theory of the classical groups can be understood natu- 
rally after continuation to negative dimensions. For example, the dimension formulas for 
SO(2n) and Sp(2m) can be related by transposition of the Young tableaux for m = — n. 
Similar relations can also be established for the Casimir operators upon interchanging 
symmetrization and antisymmetrization. These observations can be written in a compact 
form e.g. SO(—2n) = Sp(2n) where the overbar means interchange of symmetrization and 
antisymmetrization. For the Kac-Moody algebras one can now establish a relation for the 
levels by equating the central charges of the Sugawara energy-momentum tensors. For the 
case we are interested in one finds [28] : 

so(^2k) K = sp(2k)_K ■ (6.6.6) 
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Thus, from the coset realization of W£> m one finds that 

sp(2k) K <g> sp(2k)_i 

WV- k = ^ 2 - = W(2, 4, . . . , 2k(k + 2)) (6.6.7) 

sp(2k) R _i 

with a first generic null field at conformal dimension 2k 2 + 4k + 5. This fields content can 
e.g. be inferred from a character argument [28]. The algebras WE>_fc are also supposed to 
be unifying W-algebras for the WC n Casimir algebras [27]: 

WV- k ^WC n at c Cn (n + k + l,2n+2k + l). (6.6.8) 

The field content of (6.6.7) is consistent with the truncations of table 8 predicted for the 
minimal models of (6.6.8). 

To summarize, the space of all W-algebras gives rise to complicated structures in- 
cluding a 'unifying structure' that we have discussed in this section. However, there are 
indications [28] that the only rational models of unifying W-algebras are located at in- 
tersection points with the Casimir W-algebras. This implies that unifying W-algebras do 
probably not give rise to new RCFTs, and that the classification problem of RCFTs (which 
is far from being solved) might in fact be simpler than a classification of W-algebras. One 
might hope that (super-symmetric) quantized DS exhausts all possible RCFTs (with ex- 
ceptions at (half-) integer effective central charge). Thus, it would be desirable to have at 
least a good understanding of representations of non-principal DS W-algebras where very 
little is known so far. Nevertheless, the unifying W-algebras could still be of some use 
because on the one hand they might lead to conceptual simplifications and on the other 
hand not all physical phenomena are necessarily described by RCFTs. 
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... In these colleges, the professors contrive new rules and methods of agriculture and building, and 
new instruments and tools for all trades and manufactures, whereby, as they undertake, one man shall 
do the work of ten; a palace may be built in a week, of materials so durable as to last forever without 
repairing. All the fruits of the earth shall come to maturity at whatever season we think fit to choose, 
and increase an hundred-fold more than they do at present; with innumerable other happy proposals. 
The only inconvenience is, that none of these projects are yet brought to perfection; ... 

Jonathan Swift in "A Voyage to Balnibarbi" 

7. Conclusion 

In this thesis we have discussed several aspects of quantum spin models and extended 
conformal algebras. 

We have investigated an algebraic condition that guarantees integrability of quantum 
spin models - the so-called 'Dolan-Grady-condition' (also called superintegrability) . In 
one dimension we achieved a classification of superintegrable nearest neighbour interaction 
quantum spin models based on si (2). For two dimensions the only result we obtained so 
far is a non-existence statement of certain superintegrable spin-1/2 models. Nevertheless, 
this approach seems to be feasible also in higher dimensions and deserves further attention. 

In a large part of this thesis we have studied the Z n -chiral Potts model without 
necessarily demanding integrability. We have presented an argument using perturbation 
theory that shows that the massive high-temperature phases of all Z n -spin quantum chains 
exhibit quasiparticle spectra with n — 1 fundamental particles. Since the argument relies 
on perturbation theory it applies rigorously only to very high temperatures. Due to the 
perturbative nature of the details we were not able to give it any predictive power for those 
case where some of the fundamental particles cross with scattering states. For these cases 
one needs different methods, e.g. Bethe ansatz techniques [41, 103] or numerical methods 
[69] . Nevertheless, the basic idea of approximating a multi-particle state by single-particle 
states sitting on subparts of the chain might be applicable in the entire massive high- 
temperature phase. 

Our derivation of the quasiparticle picture involving n — 1 fundamental particles is valid 
also for the scaling region near the conformal point A=1,0 = <^ = O, the only non-rigorous 
part of the proof being the radius of convergence. 

Having derived such a quasiparticle picture the main open problem to obtain the scattering 
matrix of the corresponding massive field theory. 

Furthermore, we discussed a general duality property stating equality of spectra in the 
low- and high-temperature phase. Thus, the quasiparticle interpretation for the high- 
temperature phase of the general Z n -chiral Potts quantum chain can be pulled over to the 
low-temperature phase. However, charge Q and boundary conditions R are interchanged 
by the duality transformation. In particular, for periodic boundary conditions one sees 
only energy levels above the ground state that correspond to composite particle states. 

We also studied correlation functions using a perturbation expansion for the ground 
state of the Z 3 -model and numerical techniques. Although these approaches are limited 
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to short ranges, we were not only able to estimate correlation lengths in the massive high- 
temperature phase but it also turned out that the correlation functions have oscillatory 
contributions. For very high temperatures the oscillation length is proportional to the 
inverse of one of the chiral angles L ~ </> -1 - We further observed that the oscillation 
length is closely related to the minimum of the dispersion relations for general values of 
the parameters. The relation LP m ; n = 2n is valid on the lattice with a much better 
accuracy than the well-known relation £ ~ m _1 . For special values of the parameters we 
were able to derive the relation LP m [ n = 2ix from a form factor decomposition but one 
should certainly try to gain more insight into the correlation functions, in particular into 
the long-distance behaviour of the amplitude. 

We have then turned to extended conformal algebras. First, we have shown how 
automorphisms of W-algebras can be used to enlarge minimal models by imposing twisted 
boundary conditions on the chiral fields. Some of these rational conformal fields theories 
are relevant for the Z n -spin quantum chains. Furthermore, there is a striking one-to-one 
correspondence of boundary conditions for the chiral fields of a W-algebra and the quantum 
spin chain. 

Classical W-algebras are an ideal testing ground for the study of construction principles. 
Surprisingly, certain basic properties can not as easily be derived for them as for quantum 
W-algebras. We have argued that the space of fields of a classical W-algebra should 
under very general conditions decompose into su(l, 1) highest weight representations, i.e. 
it should admit a basis of quasi-primary fields and derivatives thereof. We have also 
explained why we expect that it should be possible to choose the generators primary for a 
classical W-algebra if, in addition, the conformal dimensions of all generators are strictly 
positive and the Virasoro centre c is non-zero. 

Classical cosets and orbifolds can be treated systematically. For orbifolds with respect to a 
Z2 automorphism the generators and relations are completely under control. Concerning 
classical cosets one can in general reduce the problem of finding the generators and relations 
to a problem in the invariant theory of finite-dimensional Lie algebra. Where the solution 
to the problem in invariant theory is known (e.g. for si (2) - see [124]) the space of fields is 
also completely under control. It follows from general considerations that the space of fields 
of classical cosets and orbifolds is always an infinitely generated differential polynomial ring 
subject to infinitely many relations. 

The situation is contrary for quantum W-algebras. It seems that on the quantum level 
cosets and orbifolds are always finitely generated. This difference is due to a cancellation 
mechanism between classical generators and relations. Note however that this mechanism 
has been rigorously checked only for a few cases [31, 28]. Nevertheless, there had been some 
puzzles concerning the space of W-algebras that could be explained once this difference 
was recognized. 

A further difference between classical and quantum W-algebras is that quantum W- 
algebras can truncate at certain values of the central charge because some of the generators 
become null fields. One can argue by inspection of the Kac determinant that such trun- 
cations are in fact a very general phenomenon for quantum Casimir W-algebras. These 
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truncations give rise to new 'unifying' quantum W-algebras that interpolate these trunca- 
tions of Casimir W-algebras. The notion of 'unifying W-algebras' is a generalization and 
more precise version of level-rank-duality. Unifying W-algebras can frequently be realized 
in terms of cosets, not only but also if this is expected from level-rank-duality. In the coset 
realization of unifying W-algebras non-principal DS reductions and also negative dimen- 
sional groups play a certain role. In particular, algebras of type WD-jt arise as unifying 
W-algebras for the Casimir algebras WC n . 
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Appendix A: Perturbation expansions for two-particle states 

In this appendix we explicitly diagonalize the first order of the perturbation expansion for 
two identical particles. 

Let us first consider Qi ^ Q 2 . In this case we specialize to n = 3 such that we can choose 
Qi = 1, Q2 = 2. Now, the potential V acts in the space (4.2.9) as: 

*r(Y) lit}' 2 )) P = - A (e-«*-*> 11^)) P + 2 cos (f )e"«* + *> ||4< 2 » p) 

*r(V) \\t)>*)) P = - A ( 2 cos (f )e«*+*> H^)) P + 2 cos (f ) e "*(* +f > ||^)) P ) 

1 < j < AT — 1 

9»W H^,)) p = - ^ (2 cos (f )e*(#+* > \\t]f_ 2 )) p + e«*-*> ||^ 2 » P ) 

(A.l) 

where (/ is the projector onto the space (4.2.9). Although it is not difficult to diagonalize 
(A.l) numerically for comparably long chains (e.g. A" = 100), we did not succeed in 
obtaining a closed expression for the eigenvalues or eigenvectors. 

In the second case, i.e. Q := Qi = Q2 introduce the abbreviation W by: 

-T^c cos (| - P mm , ) W\\tf Q )) p := qr{V)\\tf Q )) P . (A.2) 



In the case of two identical excitations we will also have to distinguish between even and 
odd momenta in terms of lattice sites which is conveniently encoded in the abbreviation 
dp defined by (4.2.10). The action of the potential V now is 

W\\tf Q )) P = (e-^\\tfQ)) P ) 

W\\tf Q )) P = (e^ ||f p + e-'T P ) 1 < j < [f ] - 1 

f (e 4 * \\t Q N Q 5 )) p + e-*T ||#« )) p ) , if AT odd; 

[f Hl { (e«* ||f|.« » p + ^^ e -*f ||^» p) , if iV even 

„,„ q q u f (e^||^3»p- (-l)^ll^ ))p), ifiVodd; 
If p = / v —2— —2— J 

W¥ m ))p {s^Ht^r, if iV even. 

At first sight (A. 3) looks much more complicated than (A.l). This is however misleading 
and the matrix W can be diagonalized explicitly. In order to do so, we exploit a connection 
to graph theory (see e.g. [80]). Before performing an explicit diagonalization we first present 
a graphical representation of (A.l) and (A. 3). Each vector H^ 1 '^ 2 )) p will be symbolized 
as a '•' with the index written above. The action of the potential V is symbolized by lines, 
with the square of the matrix elements (up to an isomorphism to be presented below) 
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attached to them. Assume first that we could distinguish the two flips we make. Then the 
graphical representation for the action of the potential V (or W) would be 



2 



N-2 N-l 



= (A N -i) ■ 



(A.4) 



Here '(£&)' denotes the incidence matrix derived from the Cartan matrix of a Lie algebra 
Ck- However, the states p and Wt^'Sj)) p are proportional to each other and must 

therefore be identified. Furthermore, it turns out that for N even and odd ||4' Q )) p = 

2 

vanishes identically. This already splits the graph (A.4) into two disjoint parts. Therefore, 
a graphical representation of (A. 3) is given by: 



W 



W 



2 



AT o 



■-1 



N-3 
2 



= (An 



NP 

for N even, — — odd 



W = •- 



N_ 

2 



N_ 

2 



N_ 

2 



(%) 



for N odd 



for N even, 



2tt 



NP 



(A.5) 



even. 



Fortunately, all the graphs (A.5) have norm less or equal to 2 8 ). Because the eigenvalues 
of such graphs are classified [80] we can derive the first order explicitly. 

In the case of (A.l) the situation is different. In order to simplify the discussion 
consider the case P = <j) = 0. Then one can represent (A.l) as 



V 



N-l 
■ •- 



N-2 N-l 



(A.6) 



Note that instead of drawing a closed diagram we have represented part of it twice. It is 
easy to see that the norm of (A.6) is larger than 3 (it tends to 4 for N — > oo). The absence 
of explicit expressions for the eigenvalues of such graphs prevented us from deriving an 
explicit expression for the first order of two-particle states in the Q = sector even for the 
Z3-chain. 

Let us now explicitly diagonalize the matrix (A. 3). Although this can be done for all 
four different cases simultaneously it is more illustrative to treat first one case separately. 
The simplest case is actually N even and odd. First, we perform the following diagonal 
change of bases 



B\\t) 



Q,Q\ 



')) p: =e-^-%f Q )>p. (A.7) 

If we also identify the p with the standard basis of C~ _1 the -y — 1 times -y — 1 

matrix W satisfies 



W = B 



f° 


i 





" °^ 


1 





1 







1 




'•• 








1 


Vo 







1 0/ 



(A.8) 



*) (T k ) is the Tadpole graph. 
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The eigenvalues and eigenvectors of the matrix on the r.h.s. of (A. 8) are well-known in the 
literature (see e.g. [80], example 1.2.5). We can use them to write down immediately the 
eigenvalues and eigenvectors of the matrix W: 



\r^ Q )) P : = B 



sin 



Vsin(^,)y 
W||rf'0»p = 2oofl(^)||r^»p > 



(A.9) 



N 

1< k < 1. 

~ ~ 2 



Putting things together one obtains the first order expansions (4.2.11) and (4.2.12) for 
4< N even, ^ odd. 

The same argument can be applied to the remaining cases. The only additional consid- 
eration is that we need a second change of bases M in C^ -1 . For N odd it is defined 
by 



M := 



V2 



( 1 







V-i 







1 

-1 



1\ 

• •' 



'•• 

1/ 



N 



1 x N — 1, N odd 



(A.lOa) 



whereas for even it is defined as follows 



/ 1 





V2 



1\ 

• •' 



1 





1 

y/2 
1 





V-i 



o 



'•• 

1/ 



N -lx N -1,N even. (A.lOb) 



With the definition (A. 7) of the N — 1 x N — 1 matrix B we can write W as 

••• OX 



w even 







v^ odd 



= BM 



/° 1 

1 

1 

Vo ••• 





1 





1 





1 

0/ 



M~ x B~ l 



(A.ll) 
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where W even is the [A] x [A] matrix W for even lattice momenta (5p = 1) and W odd is the 
[^T^] x [^T^] matrix for odd lattice momenta (with reversed order of basis vectors). 

Again, we can use the well-known results [80] to write down the eigenvalues and eigenvec- 
tors of the matrix W: 



( 



\t^ Q ))p: = B 



sin 



N 



\ 



sin 



\ N J 



sin ^ 



sin 



W\\t®' Q }} p = 2 cos 



sin 



(2fc-CK[f 

Wu JV 



TV 



TV 



-4= for A even 
for A odd / 



(2k -5p)tt 



A " |Tfc 



1< k < 



A 

~2 



(A.12) 



Inserting definitions one immediately obtains the final results (4.2.11) and (4.2.12) from 
(A.12). 

Let us now turn to the second order. For simplicity we restrict to the Q = 2-sector 
for n = 3 and odd A. We abbreviate the resolvent by g and its values by: 

#3,1 := -4V3cos (|) , i?i,2 := ^v^cos (^p) , Ro.i ■= W3cos (^p) • 

(A.13) 

Then we obtain for the matrix elements between the states \\t k ' )) p of the second order 
expression: 



B~ l qr(V)gr(V)qB = 



cos P 



+ 



COS T + -§ 



Ri. 



/° 1 
1 o 

o 1 



+ 



+ 



A-4 



R 



3,1 



a 



(p-f)+i 



#0,1 



1 1 

+ 



#1,2 #3,1 




1 



'•• 

1 

1 0/ 



cos 



+ 



#3, 



(1 





\0 

/0 



1 

Vo ••• 



1 







o/ 



1 

1 0/ 



(A.14) 



For a chain of length A (A odd), the matrix (A.14) is of size [A] x [A], it is straightforward 
to take the matrix elements of this matrix between the first order eigenstates (A.12). We 
omit the explicit form of this second order expression because it would be almost as long 
as (A.14). 



81 



Appendix B: Proof of the symmetries of the Hamiltonian 

This appendix contains a proof of the last line of (5.4.7). 

In order to show the third line of (5.4.7) we look at the action of the Hamiltonian (3.1.1) 
on the states (3.1.11): 

(N n-1 \ 
j=l k=l ) 

-A E "t atVte^e^tt^e*'-"-"*®^ ll«i • - - .»» , 

{ai,...,ajv} k=l 

(B.l) 

where the ||ai . . . a at)) p form a suitable basis of the space with momentum P and charge 
Q. The /fe{Vi"'°a2} an d #fe°{Vi are certain integers depending on the choice of basis. 

The yj ai '---' aN f are cer tain rea/ constants. The crucial step in the derivation of the above 
identity is to show that the basis can be chosen such that 

j.{ ai ,...,a N } = kQ -l (B , 

J k,{t 1 ,...,a N } V > 

if Q~ x G Z n . The existence of a basis ||ai...a/v)) p satisfying (B.l) and (B.2) can be 
seen as follows. Starting from an arbitrary state with charge Q and momentum P, e.g. 
the state \\s®)) p defined in (4.2.2), one can iteratively define the phase factors of the new 
basis vectors ||ai . . . ajv)) p appearing in (B.l) such that (B.2) holds. It remains to check 
that this definition is consistent which is equivalent to showing that a basis vector which 
is already introduced turns up with the proper phase factor in (B.l). Since all states 
||ai . . . ajv)) p and ||6i . . . 6jv)) p with = Xlj=i a j ~ bj have an equivalent behaviour under 
the action (B.l) of the Hamiltonian, one can restrict to the states \\s^)) p when checking 
consistency. 

Let us now consider the states \\s^)) p more closely. First we observe that one application 
of the potential V has contributions also to \\s®)) p. The two relevant terms in V||s^)) p 
are aQe lP \\s^)) p and a n -Qe~ tP \\s®)) p. These phase factors are compatible with (B.l) 
and (B.2) because e l2Pm Q is an nth root of unity. Now it remains to check that iterations 
of (B.2) that do not lead to states \\s®)) p in intermediate steps yield the correct phase 
factors. It is easy to see that the only non-trivial conditions come from iterations of type 



\\s Q ))p - ||(Q + fci)(-fci)o...))p UQ + h + k 2 )(-k 1 -k 2 )o...)) P ... » ||s«» P 

(B-3) 

without a final shift in the state ||s^)) p. The choice (B.2) for Q 1 an integer is obviously 
consistent with eq. (B.l) since X)s=i has to be a multiple of n if \\s®)) p is required to 
reappear without a shift in (B.3). This concludes the proof that one can indeed choose a 
basis such that (B.2) is valid. 
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The isomorphism stated in (5.4.7) can now be read off from the following actions of 
and (!#>)+: 




H^Wh .. .i n ))p^ q+p = - | Y.Y. akU ^ I ■■**)) p~,q+p 

{a 1 ,...,a N } k—1 

(N n-1 \ 
3 = 1 fc=l / 
n— 1 , , 

{ai,...,ajv} fc=l 

(B.4) 

In (B.4) we have used a* k = a n -k and a* k = e 2lT ' lzk ak- We would like to conclude the 
proof by noting that for very particular values of the momentum P and the chain length 
N some of the eigenstates ||ai . . . ajv)) p in (B.4) might turn out to be identically zero. In 
this case, the last line of (5.4.7) is strictly true only after inserting projection operators. 
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Translations of citations 



Aristotle, "The Physics" : 

In all sciences that are concerned with principles or causes or elements, it is acquaintance 
with these that constitutes knowledge or understanding. For we conceive ourselves to 
know about a thing when we are acquainted with its ultimate causes and first principles, 
and have got down to its elements. Obviously, then, in the study of Nature too, our first 
object must be to establish principles. 

Immanuel Kant, Preface to the second edition of "The Critique of Pure Reason": 

Mathematics and physics are the two theoretical sciences which have to determine their 
objects a priori. The former is purely a priori, the latter is partially so, but is also 
dependent on other sources of cognition. 

Werner Heisenberg: 

So far we do not know yet in which language we can talk about events in the atom. We 
have a mathematical language, i.e. a mathematical scheme, which allows us to compute 
the stationary states of the atom or the probabilities of a transition from one state to 
another. But we do not know yet -at least not in general- how this language is connected 
with the usual language. 

Niels Bohr: 

On the other hand we encounter difficulties which lie so deep that we do not have any 
idea of the way of their solution; it is my personal opinion that these difficulties are of 
such a nature that they hardly allow us to hope that we shall be able, inside the world of 
the atoms, to carry through a description in space and time of the kind that corresponds 
to our ordinary images. 

H.A. Lorentz: 

All this has great beauty and extreme importance, but unfortunately we do not really 
understand it. We do not understand Planck's hypothesis concerning oscillators nor the 
exclusion of the nonstationary orbits, and we do not understand, how, after all, according 
to Bohr's theory, light is produced. There can be no doubt, a mechanics of quanta, a 
mechanics of discontinuities, has still to be made. 
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